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ABSTRACT 



Parameter plane techniques of analysis and design are applied to 
variable parameter linear electric networks. Parameter plane theory is 
extended in several areas. 

A theorem permits the constant zeta and omega loci to be obtained in 
non-parametric form. 

The class of network functions to which parameter plane techniques 
may be applied is broadened by permitting parameters to appear non- 
linearly in the polynomial coefficients. 

Algebraic design methods are presented which permit the solution for 
parameter values which will simultaneously meet many different types of 
specifications . 

Pole-zero sensitivity formulas are derived from the non-parametric 
equations of the constant zeta and omega loci. 

Equations for curves of constant bandwidth and Q are derived which 
permit these characteristics to be displayed on the parameter plane. 

Examples are presented which demonstrate the application of these 
techniques to several common networks. 
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1 . 



Introduction 



Linear active and passive electric networks frequently have one or 
more variable parameters. In certain instances network elements are in- 
tentionally varied so as to modify the characteristics of the network. 

In other instances the network has no intentionally variable parameters, 
but due to component tolerances, aging, heating etc., variations do oc- 
cur. If it is assumed that these variations occur sufficiently slowly 
so that the network may be considered time invariant, and if it is de- 
sired to predict or control the effects of such variations on the charac- 
teristics of the network, then the methods to be described can be applied. 

Parameter plane techniques were first introduced as a technique for 
analysis and synthesis of feedback control systems. They were used to 
investigate the effects upon the roots of the control system character- 
istic equation as two parameters within the system were varied. 

In the early 1950*3 Mitrovic (1) considered system parameter varia- 
tions which directly affected two of the coefficients of the system 
characteristic equation. For example: 

a^s + a^ ^s 11 ^ + ...+ a^s^ a 2 S ^ + ^l s * = ^ (1^4) 

where: a^ are constant and real 

Bq and are variable and real 

Mitrovic* s original work resulted in curves of constant zeta, omega, and 
sigma which he drew upon the B^-B^ plane and which specified the roots of 
Equation (1-1) for any choice of B^ and B 1 . 

In 1964, Siljak (2) extended Mitrovic* s work to include cases where 
variable system parameters, say o< and (3 appeared linearly in one or 



more of the coefficients of the characteristic equation. Specifically 
Siljak investigated the case where: 



where: 



£ 

I - A 



V 



= o 



( 1 - 2 ) 



a . = b ,°< +c.B + d. 
11 l ' l 



and: 

\ 

b . , b., d. are constant and real 

l l l 

c* , (3 are variable and real 

Siljak, like Mitrovic, obtained curves which were plotted upon the 
plane and which specified the roots of Equation (1-2) for any choice of 
and (3 . It was Siljak who first referred to such graphs as "parameter 

planes". 

The parameter plane provides information about the roots of the 
characteristic equation similar to that provided by root loci, but unlike 
root locus techniques, which investigate the effects of varying one system 
parameter, parameter plane techniques are inherently best suited for in- 
vestigating the effects of varying two system parameters. 
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2 . 



Basic Parameter Plane Techniques 



Parameter Plane techniques evolved from a study of the characteristic 
equation of control systems and of how the roots of such equations were 
affected as two parameters of the control system were simultaneously varied. 
More basically however, the parameter plane is a method for investigating 
the roots of variable coefficient polynomials, and hence is applicable in 
a much wider sphere than that of control systems alone. This chapter will 
review the parameter plane techniques as presented by Siljak (2) and will 
point out some of the reasons that have limited the scope of applicability 
of these techniques. Subsequent chapters will present new work which 
overcomes some of these limitations. 

Let us consider a polynomial equation: 

n V 

Hs) = Y1 V = 0 <2-D 

k = o 

where the coefficients a k (k = 1> . . . ,n) are real and s is the 

complex frequency variable s = cr + jtO = - 1 -p z 

If one forms powers of s it can be shown (2) that the following relation 
holds: 

s k = co k ( T k (-y ) + j / 1- U k (- } )) (2-2) 

where: 

V f > = ( - i)k v ? > 

(2-3) 

y- y ) = (-D k+1 y y ) 

The U k and T^ are Chebyshev functions of the first and second kind 
respectively and are given by the following recursion relations (2): 
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T k + 1<? > - 2 ? V J > + Vi< ? > ■ 0 

(2-4) 

"k + 1< y > - 2 ? v ? > + \.i< ? > - ° 

where: 

V f ) - !->. V f > - ? . V ? > * 0 ’ B l ( ? > * 1 



Tables of these Chebyshev functions are included in Appendix III. 

Substituting Equation (2-2) into Equation (2-1) and setting the real 
and imaginary parts of the resulting equation equal to zero independently we 
obtain: 



1 

kro 

£ 

K=0 



a k U.\( - ?> - 0 

a k V -f > - 0 



(2-5) 

(2-6) 



Employing Equations (2-3) and (2-4) we obtain these expressions in 



terms of the tabulated Chebyshev functions: 

n 



Z<- i)ka k"Vi ( ? > - ° 

K-O 

£ 

K=0 



(-l) k a k 6) k U k ( £ ) = 0 



(2-7) 

(2-8) 



Let us consider, as Siljak did, that the coefficients a k of Equation 
(2-1) are linear functions of the variable parameters o< and (2> as follows: 



where: 



\ = V + c kP + d k 



(2-9) 



b^, c k , d^ (k = 0, 1, ...,n) are constant and real 

. i 

o( , (3 are variable and real 
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Substituting Equation (2-9) into Equations (2-7) and (2-8) we obtain: 



B 



+c iP 



V< +C 2 



P 



where: 



b i - 



'■AA-i 

i 

C l= £(- 1 ) k c k w\.i 
KrO 
n 



<=<> 



- - D, 



- D, 



( 2 - 10 ) 



B 2 = 



C 2 = 



°2 = 



2 <-»vn 

K-O 

2 < ; l) ’V' u i 

K=0 

t- 



dVA 

k k 



( 2 - 11 ) 



Equations (2-10) are two linear equations in the variables °( and ^ . 
Solving these equations we obtain: 



* 



C 1 P 2~ C 2 P 1 

B 1 C 2- B 2 C 1 



P 



B 2 D r B i D 2 

B 1 C 2- B 2 C 1 



( 2 - 12 ) 



Equations (2-12) give o< and in terms of ^ and *0 . Thus if we 
fix ^ and vary k) over the range 0 < 00 we will ob- 

tain a curve in the o( - (S plane which specifies the relation between o( 
and necessary to cause a pair of complex roots of Equation (2-1) to 

have a damping ratio ^ equal to ^ Curves resulting from fixed 

values of ^ are referred to as constant zeta loci. If we fix ^0^ 

and vary ^ over the range -1 ^ ^ K ^1 we will obtain another 

curve in the o( - j3 plane which specifies the relation between o( and ^3 
necessary to cause a pair of complex roots of Equation (2-1) to have an 
undamped natural frequency U) = u) g. Curves drawn for fixed values of 
k) are referred to as constant omega curves. 
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If we substitute s = CT in Equation (2-1), where sigma refers to 
values of s along the real axis in the s-plane, then we obtain: 



n 

2 a k cr k = 0 (2-13) 

K=0 

If we now substitute for the linear relation (2-9) we obtain: 

b k °' k + i 8 i c k°' 1 ' + ^ d k< rk -° < 2 - i4 > 

K-0 K-0 KrO 

For fixed values of sigma J Equation (2-13) represents a straight line in 
the c< - plane. This straight line specifies the relation between <*< 
and |3 necessary for Equation (2-1) to have a real root at s = <T . Such 
lines drawn upon the c< - p plane are referred to as constant sigma loci. 
If the reader is interested in the mapping which occurs between the 
s-plane and the parameter plane or in the properties of this mapping then 
the author refers him to reference (2) where this is discussed in detail. 



Example I : 

Let us consider that we have constructed a parameter plane for some 
third order polynomial and that this appears as shown in Figure 1. 

Let us consider that parameter values o< and |3 are chosen to be c< 
and (3 q respectively, and that at these coordinates on the parameter 
plane it happens that the \ anc * ^2 '*' OC ^ intersect * We may 

thus assert that for this choice of parameter values the polynomial has 
the roots: 



s 



s 



<r 



1 




j <0 2 <1- ?/>* 



Let us now presume that we wish to change the damping ratio to 



r 



o* 
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- 




Fig. 1. Example of a Parameter for some 3rd Order Polynomial 

but wish to leave the undamped natural frequency unchanged at ^ 2 * 

We therefore find the intersection of the ^ ^ and U ) , loci. In this 
example this occurs at = o( (3 = /3 Adjusting the para- 

meters and f3 to these values will guarantee complex roots at s « 

2 ^ 

- J Q u) 2 ± 2 10 2 (1- $ q )' 2 . The real root which was at s = O' ■ 

has moved to a new location between O' ^ and O' 

Thus we see that the parameter is a general tool of both analysis 
and synthesis in two parameter systems. 
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3. Scope of this Work 

This section outlines the general content of this manuscript and 
introduces new work which overcomes some of the limitations of parameter 
plane theory which have heretofore prevented its application to electric 
networks. 

General Content 

In this work the author investigates the application of parameter 

plane techniques to variable parameter electric networks. The objective 

* 

of this work is to provide means of analysis which will show how the net- 
work characteristics change as the parameters are varied and to provide 
means of design which will permit values of the adjustable parameters to 
be chosen so as to cause the network to have desired characteristics. To 
this end both graphical and analytic techniques are presented. The net- 
work characteristics which can either be determined or controlled are 
pole and zero locations, frequency response, Q and root sensitivity. 

Section 4 through 8 develope the mathematical theory which related 
the network characteristics to the variable parameters. The mathematical 
theory presented in these sections forms the basis for both the graphical 

4 

4 4 

and analytic techniques. 

Sections 9, 10, and 11 contain examples of graphical analysis and de- 
sign of selected networks which have historically been of considerable 
interest. These examples provide not only an illustration of the poten- 
tial of the parameter plane techniques, but also provide normalized results 
which will aid designers in employing these networks. 

Section 12 considers several other networks of general interest. 

These networks are considered in enough detail to show the applicability 
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of parameter plane techniques to their analysis and design. 

Section 13 presents an algebraic method of design which employs the 
parameter plane equations developed in earlier sections. These algebraic 
methods of design can be applied to any problem which may be graphically 

designed and additionally to the much larger class of problem which have 

/ 

more than two parameters. There is theoretically no limitation on the 
number of parameters which may be simultaneously considered when using 
these algebraic techniques of design. 

Section 14 concludes the manuscript with the author 1 s comments on 
the work and with his recommendations for additional work in the area of 
network analysis and design by parameter plane techniques. 

Limitations of Earlier Parameter Plane Theory and Introduction to New Work 
Parameter plane theory as it has heretofore existed has several 
limitations which have restricted its application to electric network 
analysis and design. The remainder of this section discusses these limita- 
tions and introduces new work aimed at overcoming them. 

Form of the Coefficients 

Consider the restrictions which have been placed on the polynomial 
coefficients - namely that the variable parameters must appear linearly in 
the coefficients: 

a k ■ b k* + c k P + d k 

This restriction was not too severe when parameter plane techniques were 
applied primarily to feedback control systems since many useful systems 
have characteristic equation coefficients of this form. Electric networks 
on the other hand, seldom have coefficients of this type. Consider for 



9 



example the Parallel -T network shown in Figure 2: 



■R "R 




Fig. 2. Symmetrical Parallel-Tee Network with Variable Elements 
This network has the transfer function: 

E q (p) = o<(3p 3 + 2 /3 p 2 + 2/9p + 1 (3- 

E i (p) oC/S P 3 + (2 /3 + 2 */3 + ©< )p 2 + (2 /3 + <* + 2) +1 

where: p = RCs 

Observe that Equation (3-2) contains terms o((3 as well as the linear 
terms permitted by Equation (3-1). Thus we are unable at this point, due 
to the linear requirements placed on the coefficients, to determine by 
parameter plane techniques how the singularities of (3-2) are altered as 
the parameters U and /3 are varied. New work to be presented in 
Section 4 will overcome this difficulty. 

Constant Zeta and Omega Loci 

As derived by Siljak (2), the constant zeta and omega loci are 
presented in parametric form as functions of the independent variables 
omega and zeta respectively. This may be seen in Equation (2-12) where 
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the constant zeta loci are obtained by fixing zeta at selected values and 
generating the loci by varying the independent variable omega over the 
range 0 to infinity. Similarily, the constant omega loci are obtained 
by fixing omega at selected values and then varying zeta as an independent 
variable over the range -1 to +1. 

Clearly it is desirable to obtain the equations of these loci in non- 
parametric form. By achieving this we simplify the plotting of these loci, 
and more important we obtain these loci as families of plane algebraic 
curves in the o< - plane. This will permit the mathematical procedures 
for the investigation of plane algebraic curves to be applied when analyz- 
ing the behavior of these loci. This is an important step if procedures 
for sketching the constant zeta and omega loci are to be developed. In 
Section 5, procedures are presented for reducing these equations to non- 
parametric form. 

Frequency Response Characteristics 

Since frequency response characteristics are often used in specifying 
the behavior of an electric network it is essential that any method which 
purports to be a general tool for analysis and design of such a network 
include procedures for determining and controlling the essential features 
of its frequency response. Choe (3) developed formulas for displaying 
curves of constant bandwidth upon the Mitrovic plane, but these techniques 
were not applicable to the parameter plane. In Section 6 formulas are 
developed for displaying curves of constant bandwidth upon the parameter 
plane. These formulas are not only applicable to electric networks, but 
to any rational system function which can be displayed upon the parameter 
plane. 
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Constant Q Loci 



Second order network transfer functions can be characterized by their 
Q, much in the same way that the characteristic equation of a second order 
servomechanism is characterized by its damping ratio. Procedures are 
developed in Section 7 for displaying loci of constant Q upon the parameter 
plane. These specify the value of Q which results for any choice of network 
parameters or alternately specify in what fashion the network parameters 
must be changed in order to obtain some desired value of Q. 

Root Sensitivity to Parameter Variations 

The sensitivity of the roots of a polynomial to changes in parameters 
contained in its coefficients can be catagorized as macroscopic or microscopic. 
Macroscopic sensitivity refers to the sensitivity of the roots to large changes 
in the parameters. Microscopic sensitivity refers to infinitesmal changes in 
the parameters. The constant zeta and omega loci drawn upon the parameter 
plane permit macroscopic sensitivity to be determined by inspection. Micro- 
scopic sensitivity may be determined from formulas derived by Kokotovic and 
Siljak (5) from the parameter plane equations or may be computed directly 
from the non-parametric equations of the constant zeta and omega loci which 
are developed in Section 5. Section 8 describes how both macroscopic and 
microscopic root sensitivity to parameter variations may be determined. 

Examples of Parameter Plane Techniques Applied to Networks 

Sections 9-11 contain detailed examples of parameter plane techniques 
applied to specific networks. These examples were carefully chosen in 
order to both fully demonstrate the current state of parameter plane tech- 
niques, and at the same time to produce results useful to persons desiring 
to employ these networks. The results which are presented are general and 
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are presented in normalized form. 

Section 12 contains examples of several networks with variable para- 
meters which were studied only so* far as to show the applicability of 

* t 

parameter plane techniques to their analysis and design. 
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4. 



Parameter Plane Equations for Coefficients of the Form: 



a k ‘ V* + °k « + \ + f k 

In this section equations will be obtained which will permit para- 
meter planes to be constructed for polynomials whose coefficients are of 
the form: 



-v* + ' k e +d k^ +£ k 



(4-1) 



where: b^ 5 c^> d^ are constant and real 

c< , (3 are variable and real 
Consider the polynomial equation: 



n 



s =0 



F(s) ) a k 

KiO 

Let the coefficients a^ be of the form specified by Equation (4-1). 



(4-2) 



Let : 



= *> k [T k (-f ) + j Jl- f 2 U k (-J)] (4-3) 



where: and are Chebyshev functions discussed in Section 2, 

Substituting (4-3) into (4-2) we obtain: 
n 



n i n n 

F ( s) . x a k T k ( - ? > ± j 7 i - r 2 2 a k , ° k v-f > 

K=0 K=6 (4-4) 

Equating the real and imaginary parts of Equation (4-4) to zero independent- 
ly we obtain: 

n 



z \« kT k<-? > ■ ° 



K=o 



V-p-0 



(4-5) 

(4-6) 



K-O 



14 



but from Equations (2-3) and (2-4): 



V? > ■ <■') Vf> 



\<- ? > - <-l )k+l \< f ) (4-7) 

"k-^- ? > - <‘ 1 > k "k-^ f > 

Substituting Equations (4-7) into Equations (4-5) and (4-6) we obtain: 



(-l) k a «J k T < ? ) » 0 



(4-8) 



k^O 



2 <- i >\“ k v? >-° 

K-C 



(4-9) 



Substituting for the coefficients in Equations (4-8) and (4-9) the 
form specified by (4-1) we obtain: 

* x <-»\ i 3 y <- i ) k -=k “V f > + 

K=0 K=0 

«|9 ^ <- l >\ " V ? > + £ <- l > kf k "V f > - 0 

K = d AC=P 

(-1)\ to\(f ) + /3^ (-1)\ w k u k (f ) + 

(-1)\ AJ k U k ( f ) + ^ (-l) k f k uJ k U k ( £ ) - 0 



(4-10) 



K-0 



which may be written: 



B x o< + C x /8 + D x o< (3 + F - 0 

B 2 o< + C 2 0 + D 2 »(/(3 + F 2 = 0 



(4-11) 

(4-12) 

(4-13) 



where: 



°i- 



<- i >\ " k v r > 



n 



c i = i ( - i) \ ? > 



b 2 - 



C 2 = 



X < - i)kb k« k v? ) 






(-l) k C 4J k U k ( f ) 



*=o 






(4-14) 
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D 1 ■ £ <- 1 >\ "\< V °2 ■ 2 ( - 1)kd k > 

K^o >r«o 

F 1 ■ $ <_1)k£ k "V ? > f 2 * $ <- 1)kf k "\< i > 

KiO KaO 

Eliminating o< from Equations (4-12) and (4-13) we obtain: 



F 1 + C ifl 
B ! + D 1 /6 



F 2 + C 2 £ 



b 2 + ° 2 /3 



= ©< 



(4-15) 



Which may be written: 



AcD ^ 4 ^App 4 A eg) ^ + Aptf ° 

r-8 



(4-16) 



Where : 



CD 



C 1 °2 



D 1 °2 



^FD 



F F 

1 2 



D 1 °2 



, etc . 



Equation (4-16) has the solutions: 

"(ApD + A CB ) — ^ + A CB ) ^^CD^FD 



P 



(4-17) 



2A 



CD 



>9 



Let us define: 



R p “ ^Ap D + A CB ) ' ^Aqd^FD 

This permits us to write Equation (4-17) as: 



(4-18) 



P = 



A 



CD 



[ -<A FD + A™> ±R/8 1 



CB' 






(4-19) 



Proceeding in a similar fashion we may eliminate ^3 from Equations (4-12) 
and (4-13) to obtain: 



r t + b i0 ( 

c i 4 m 



F 2 4 B 2 °< 
C 2 4 D 2°< 



As before, we may rewrite this as 
. 2 



(4-20) 



A 



ED 



4 <A B c 4 A F D )o < 4 A 



FD 



FC 



(4-21) 
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we obtain: 



Solving Equation (4-21) for 

« - 2^7 [ -<A pd + A bc > ±R* ) W-22) 

Where we have defined: 

R o< = + ^FD^ ” ^^BD^FC (^ -23 ) 

2 2 

We will now show that = R^g 

From Equation (4-23) we write: 



R^= (B 1 C 2 -B 2 C 1 + D 2 F 1 - D 1 F 2 ) 2 - 4(F 1 C 2 - F 2 C 1 )(B 1 D 2 - D^) (4-24) 

Expanding Equation (4-24) we obtain: 

R « ■ B 1 C 2 + B 2 C 1 + B 2 F 1 + B 1 F 2 + 4 B 1 C 1 D 2 F 2 + 4 B 2 C 2 D 1 F 1 ' 2 B 1 B 2 C 1 C 2 ' 

2 B 1 C 2 D 2 F 1 - 2 B 1 C 2 D 1 F 2 - 2 B 2 C 1 D 2 F 1 ‘ 2 B 2 C 1 D 1 F 2 ‘ 



2 D D F F 
12 12 



(4-25) 



Expanding Equation (4-18) for R^ we obtain the same expression as in 
the right hand side of Equation (4-25). Thus: 



R 2 - R 2 * t R 2 



(4-26) 



Thus we may state from Equations (4-19) and (4-22): 





« - 2 Abd (-<A fd *A bc > ±») 


(4-27) 




^ ~ 2 A rn ^^FD + ^ CB^ - 

(jU 


(4-28) 


Noting 


that Aiir ~ " /\n-o we may rewrite Equations (4-27) 

dL 1 UJ 


and (4-28) 


as : 






o< = 


Q 

2S; d 1 -(An,* A bc ) + M*) i *-(i 


(4-29) 


|8 = 


2 Acd ' ' <Afb + Abc> + < ' 1>mRl 


(4-30) 
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Equations (4-29) and (4-30) define four possible solution pairs depending 
on the values of k and m. 

Let us now determine which of these possible solutions are actually 
solutions to Equations (4-12) and (4-13). Multiplying Equation (4-12) by 
D 2 and Equation (4-13) by -D^ and adding we obtain: 

AbD°^ + AcD @ + ^ FD “ ° (4" 31 ) 

Equations (4-31) must be satisfied by any of the pairs (4-30) which are to 
be solutions to Equations (4-12) and (4-13). Substituting Equations (4-29) 
and (4-30) into Equation (4-31) we obtain: 

- A FD - ABC + ( -l) k R - AfD + AbC + (-l) m R + A™ = 0 
2 2 2 2 2 2 



(4-32) 

Which after simplification becomes: 





<-D k + 


(-D m = o 


(4-33) 


In order for (4-33) 


to be satisfied it is necessary that 








k ^ m 


(4-34) 


Thus we may rewrite 


Equations 


(4-29) and (4-30) as: 




°< k 


1 

2 Abd 


I -<A fd + Abc> + <-') kRI 


(4-35) 


£k 


2 A CD 


t -<A fd + A bc > + <-D k+1 R] 


(4-36) 



Where k = 0 and k = 1 generate the two solution pairs of Equations (4-12) 
and (4-13) . 

Equations (4-35) and (4-36) may be used to program a digital computer 
to generate the constant zeta and omega loci. The basic steps to be pro- 
grammed are: 
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Constant Zeta Loci 



1. Fix ^ ~ $ i* This makes Equations (4-14) functions only of 

the independent variable omega. 

2. Increment omega and for each value of omega compute the solu- 
tion pairs ( ^ ^ k = 0, 1 using Equations (4-35) and (4-36). 

3. Connect all solutions ( o( ^3^) t ^ ie or ^er of increasing 
omega. This generates one branch of the constant zeta locus for 

4. Connect all solutions ( t ^ ie orc ^ er of increasing 

omega. This generates the other branch of the constant zeta locus for 



I" ?!■ 



5. Repeat steps 1-4 for ^ * ^i+l J = ^ i+2 > etc# 



until as 



many constant zeta loci are obtained as desired. 



Constant Omega Loci 

1. Fix Cd - cd y Thus makes Equations (4-14) functions only of 
the independent variable zeta. 

2. Increment zeta and for each value of zeta compute the solution 

pairs ( ^ ^) ; k = 0, 1 using Equations (4-35) and (4-36). 

3. Connect all solutions ( °^ q> fS q) in the order of increasing 
zeta. This generates one branch of the constant omega locus for t0 = u) 

4. Connect all solutions ( o( ^) in the order of increasing 

zeta. This generates the other branch of the constant omega locus for 

U) = U) j . 

5. Repeat steps 1-4 above for 10 = b) = U) . t etc. until 

J+l J+ 2 

as many constant omega loci are obtained as desired. 
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5. 



Non-Parametric Equations for the Constant Zeta and Omega Loci. 



In this section procedures are presented which show how the equations 
of the constant zeta and omega loci, regardless of the form of the co- 
efficients a^, may be reduced to non-parametric form. 

As developed in Section 4, the parameter plane equations for co- 

•f f, take the form: 
k 

c< = A( g t Cd ) (5-i) 

£ - B( ^ , u)) 

The constant zeta and omega loci are obtained from these equations 
by fixing one of the variables, zeta or omega, and considering the other 
variable as independent. Varying this independent variable over the range 
of interest generates the desired locus. Our objective in this section 
is to obtain the equations of the constant zeta or omega loci from (5-1) 
by eliminating omega or zeta, respectively. Achieving this will produce 
the non-parametric equations of the desired loci: 

z( * , ft , $ ) = o (5-2) 

W( Of , ft , Cd ) = 0 (5-3) 

Two theorems, basic to this development, are proved in Appendices I 
and II. These are: 



efficients a^ = b^o( + c^ (3 + d^ o (jQ 



Theorem I 

The Chebyshev functions U^( 

V ? > Vi ( ? > - ”i-i< ? > 



^ ) obey the following relationship: 

VP = " u i-j ( P (5_4) 



for all integers i and j. 
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Theorem II 



If: 







Z <- 1)ka k 




and if: 




KiO 

fi 








y f > - ° 




Then: 




KzjO 








Z 


(5-5) 






K=0 




for all integers k and 


j. 






Equations of the Constant 


Zeta Loci 




Using Theorem II 


for 


the cases j = 0, j = n we obtain: 




f(oj ) 


= 


Z ( - 1)ka k > ■ 0 


(5-6) 






KzC 




g(^ ) 


= 


2 ^ Vn< r > - ° 


(5-7) 






K-0 





Equations (5-6) and (5-7) can be considered to be polynomial equations of 
order n in (0 . We wish to eliminate (/) from these equations in order to 



obtain Z(a^, ^ ) = 0, which is the non-parametric representation of the 
constant zeta loci. This may be accomplished by employing Sylvester's 
Dialytic Method of Elimination (6). The basic steps are shown below for 
two polynomial equations of differing order in tO , from which (0 is 
to be eliminated. 

Given the polynomial equations: 



f( CO ) = f co m + 

m 


f „-l “j"' 1 + ••• + £ 0 * 0 


(5-8) 


g( ti) ) = g n M) n + 


Sn.luT 1 + ... + g Q - 0 


(5-9) 



The quantity (j) may be eliminated from Equations (5-8) and (5-9) by the 
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. . . , 1 in turn. 



following sequence of operations: 

1. Multiply (5-8) by 60 ", 0 ) 

2. Multiply (5-9) by GO* \ GO™ 1 in turn. 

3. Form the detenninant of the coefficient matrix of the 

resulting set of (nfm) equations which are linear and 
homogeneous in the (n-rm) quantities 60 * 1 * “ \ UJ 1 
CO , 1. 



4. This determinant, when forced equal to zero, establishes the 
relation among the f^ and g^ necessary if Equations (5-8) 
and (5-9) are to have common roots. 

The resulting determinant has the form: 



<w - 



-T» 

'a 


• 


• • 


f 2 


X n 

-L 


-T» 

x 0 


0 


• 


• 


• • 


• 


• 


• 


0 




0 


~n 


• • 


♦ 


^2 


-T* 

"1 


f 0 


0 


• 


« • 


• 


• 


• 


0 




0 


0 


X • 

m 


• 


• 


x 2 


f l 


jy 

x 0 


0 


• • 


• 


• 


• 


0 




• 

• . 




\ 






















0 




0 
















0 


n 


• 


jy 

x 2 




wr* 

-o 


* 0 




• 


• • 


s 2 


Sl 


s o 


0 


• 


• 


• • 


• 


• 


• 


0 


(5-10) 


0 


rr 

o n 


• • 


• 


S 2 S X 


s o 


0 


• 


• • 


• 


• 


• 


0 

■* 




• 

0 


• 


• • 


• 


• 


• 


• 


0 


S n 


• • 


• 


S 2 


Si 


So 





In the case of Equations (5-6) and (5-7), the constant zeta loci can 
be obtained from Equation (5-1G) by making the following substitutions: 



£ k * <Avf> 

% ■ 



(5-11) 

(5-12) 
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The following example w ill demonstrate the technique: 



Example 2 

Consider again the symmetrical parallel-tec network shown in Figure 
2, Section 3. For the zeros of the transfer of this network. Equations 
(5-6) and (5-7) become: 

o "> 

(5-13) 



- o<p <^U 3 ( £ ) + 2 (3 <0 2 U 2 ( 3f ) - 2 f ) + U 0 ( £ ) = 0 

-0^3 cd 3 U 0 ( ? ) + 2/5 <0 2 U_ X ( ? ) - 2/BU_ 2 ( ? ) + U_ 3 ( ^ ) = 0 (5-14) 

2 ^ 

Multiplying these equations by CO , CO i an ^ 1 and forming the determinant 
Z( o( , ) , we obtain: 





-? u 3 




“^1 


fU o 


0 


0 






0 


-^u 3 


t2^2 




+ u o 


0 




z(*,/S ) = 


0 


0 


-^u 3 


v2p 2 


-2^ 


• u o 


- 0 


-• ‘Fo 


v2&U_i 


-2^U_ 2 


_ « TT 

N -3 


0 


0 




0 


“^ U 0 




-2p_2 




0 


(5-15) 




0 


0 


“^O 


• 2 ^ U _1 


- 2 ^ U _ 2 


+U 




Substituting the 


appropriate 


functions 


of zeta 


for the U, 
k 


in Equation (5-15) 



and expanding, we obtain the equation for the constant zeta loci: 



<*|6 3 (4^ 2 -l) 3 [(-32 * ^ 2 +16)/9 2 + (64* ? 4 +l6* ^ 2 - 32 £ 2 - 8oO|0 + 

(- 64 * 2 f 6 - 1 - 48®<“^ 4 - 12 o < 2 |' 2 + ex. 2 )] = 0 



(5-16) 
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When the two equations (5-8) and (5-9) are of the same order so that 
m = n, Sylvester* s (2n x 2n) determinant can be reduced to an (n x n) deter- 
minant by using Bezout*s Method of Elimination (7), To demonstrate Bezout's 
method, consider Equations (5-8) and (5-9) with m * n = 3. This pro- 
duces the equations: 



Then: 



f( CO ) = f 3 u> 3 + f 2 Oi 2 + u) 1 + f Q = 0 

g( k) ) = g 3 4) 3 + g 2 0) 2 + g 1 AJ 1 + g 0 = 0 



f 3 g - g 3 f - o 



(f 3 &) + f 2 )g - (g 3 to + g 2 )f = o 

(f 3 10 2 + f 3 (0 + fj)g - (g 3 U) 2 + g 2 U) + gj)f = o 



can be written in matrix form: 

r 

u) 

u) 

i 



(f 3 g 2 ) 


(f 3 g l ) 




( f 3 g 0 ) 


( f 3 g l) 


(( f 3 g 0 ) + ( f 2 8 




( f 2 g 0 ) 


(f 3 g 0 ) 


(f 2 g o^ 




(f i 8 o ) 


) denote 


(2 x 2) determinants 


« 

• 




( f i g j) 


f . 

i 

g i 


f . 
J 

g j 



= o 



(5-17) 

(5-18) 

(5-19) 

(5-20) 

(5-21) 



(5-22) 



(5-23) 



In order to eliminate Ol from Equations (5-17) and (5-18) by Bezout's 
method it is necessary to force the determinant of the coefficient matrix 
in Equation (5-23) to zero. Thus the constant zeta loci for the case 
m = n = 3 are given in non-parametric form by: 
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i 



(f 3 S 2 ) 


(f 3 gl> 


( f 3 g Q ) 








(( f 3 g 0 ) + (f 2 sp) 


(f 2 s 0 ) 


= 0 


(5-24) 


(f 3 g 0 ) 


( f 2 S 0^ 


( f lg 0 ) 







For cases where m = n y 3, a determinant similar to (5-24) can be 
developed by forming n equations as in (5-19) - (5-21). When m = n « 4 
such a procedure gives the following determinantal equation for the constant 
zeta loci: 



z 4 «*,/s,y> 



<V 3 > 


(f 4 6 2 ) 


(f 4 s l> 


(f 4 s 0 ) 


\ 


((^2j.S^) * ("362)) 


!(f 4 s 0 ) 5- (f^)) 


(* 3 s 0 ) 


< i 


(( f 4 S 0 ) + (f 3 S 1 )) 


((" 3 s 0 ) (f 2 s 1 >) 


{ Vo } 


;^s 0 ) 




( A 2 S 0^ 


(f i S 0> 



The following example will demonstrate how Bezout's method is applied 
to obtain the constant zeta loci: 

Example III 

Let us again consider the symmetrical parallel-tee network. Equa- 
tions (5-13) and (5-14) yield: 

f 3 = -'tfi U 3 £ 2 ' +2 0°2 f l ’ ' 2 ^l f 0 ' +U 0 
3- - C 0 s 2 = +2 (S U_ l 8l = "2(3 g 0 - +U.3 
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Substituting the appropriate functions of zeta into these expressions and 
forming the (2 x 2) determinants for the (f^g^) we obtain: 

(f 3 g 2 ) = +2<*/3 2 (4f 2 -l) 

(f 3 g 0 ) = 2 f (V 2 -d 

<f 3 g 0 ) = P? 2 - °< P )(4^ 2 -l) 

(f 2 g x ) = 4f3 2 (4^ 2 -l) 

(f 2 s 0 > - (4? 2 -ii 

(f lg() ) = +2 p (4^* 2 -l) 

When these are substituted into Equation (5-24) we obtain: 



Z 3 ( * , ft , f ) = o<(3 3 (4^ 2 -l) 3 



+2(3 -4(Sf 

-4°<p^ (4°£ £ 2 + 4(3 - * ) 

(4*^ 2 -oO -4 ? 



(4 f 2 -l) 



-V 



= 0 



+2 



Which when expanded gives the constant zeta loci in non-parametric form: 

Z 3 («(, / 8,^)» <=<(<3 3 (4f 2 -l) 3 [(-32* ? 2 +l6)/3 2 + (64* f 4 +l6 * f 2 -32 £ 2 - 

8*)/3 + (-64* 2 f 2 + 48<* 2 ^ 4 -12®< 2 Jr 2 + <* 2 ) ] = 0 
Observe that the determinant obtained by Bezout's method has the op- 
posite sign from that obtained by Sylvester's method. Sylvester's method 
produces the resultant of the functions f( £0 ) and g( £0 ), while Bezout's 
method produces the negative of the resultant (8). Since multiplying both 
sides of the determinantal equations by -1 does not introduce extraneous 
common factors of f( U) ) and g( flj ) , this sign difference can be ignored. 

It would appear at first that Bezout's method would be easier to apply 
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than Sylvester* s method because of the reduction in the order of the deter- 
minants that it affords. However, Sylvester* s determinant has many null 
elements as well as functionally simpler elements than Bezout*s determinant. 
It is the author* s opinion that Sylvester* s method is somewhat easier to 
apply. 



Equations of the Constant Omega Loci 

Using Theorem (II for the cases j = 0, j = 1, ..., j = (n-1) we 
obtain: 



n 

*Z (-dWp - ° 

K=0 

2 - 0 






(5-26) 



k^O 



These equations may be written: 



(-l) 3 a c <0°U 0 + (-l) 1 a 1 0 1 U 1 



.. + (-l) n a to n v = 0 

n n 



(-l) 0 a A <d°U + (-l) 1 a 1 «>V + ... + (-l) n a uA . = 0 (5-27) 

0-1 10 n n-1 



(-l)^a A 0)^UjT (-l)"*’a 1 O^U / • +(-l) n a^ a/ 1 U 1 - 0 



'i-f 






n '1 

Making use of the relation U(^ ) = - f ) and = 0 and rearranging 
the set of equations (5-27) we obtain: 

,1 



(-l)*a n fe3 



[(-!)% « 2 -(-l)°a io°3 ... 
2 0 

[(- l) a a « a -(-l) a " 2 a 
“ n n-2 

. . n-1 , . n-1 

-1 a “ 

v-i 

^ -t 



n v? 

(-1) a fcT 
n 



0 

0 



< -1 ) °^o^ C 





J 

t - 




U 2 




U 3 




' Tj n-i 



= o (5-28) 
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The set of equations (5-28) are linear and homogeneous in the quantities 

U, . Thus, for a non-trivial solution to exist, the determinant of the 
k 

coefficient matrix in (5-28) must vanish. This condition specifies the 
relation nocoaoary among the a„^ and *0 which apocifies the constant omega 
loci : 



W(a k ,< 0 ) = 



(-1) a^w 1 

[(-l) 2 a 2 w 2 -(-l)°W U j ... 

L(-i) n a r /-(-i) n - 2 ^:?. 

, _ v n-1 n-1 

(-D 



(-l) n a iti“ L 
n 



0 

0 

(-l)°a 0 to C 



= 0 



(5-29) 



Example IV 

Let us apply this procedure to obtain the constant omega locus for 
the zeros of the symmetrical parallel-tee network. From Equation (3-2) 
the zeros of the transfer function for this network are the roots of: 

c^|Sp 3 + 2 (9 p^ -f 2^>p 4-1 = 0 (5-30) 

Applying theorem II, we may write as in Equation (5-26): 

(-l) 1 (2^)<iu i * (-1) 2 (2 (3 ) t 0 2 U 2 + (-1) 3 (<*(3 ) = 0 

(-1)°(1) u)“U_i * (-D 2 (2/3) C0 2 U x + (-1) 3 ( <*/3) 4) 3 V 2 =0 (5-31) 

(-1)°(1) «0°U + (-1) 1 (2 + (-l) 3 (^)cJ 3 U 1 - 0 

Using the relation U , = - U, these equations can be written: 
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(-l) X (2p )tO X 


(-1) 2 (2£ )<0 2 


(-1) 3 ( )<o 3 




U 1 


[(-1) 2 (2^> ) A) 2 -(-l)°(l) u)°] 


(-i )h*p )*> 3 


0 




U 2 


t(-l) 3 (^P )o) 3 -(-l) 1 (2(3 )40 1 ] 


-(-1)°(1)40° 


0 




U 3 



Setting the determinant of the coefficient matrix to zero we obtain: 






-2$W +2(3 

( 2 {$ oJ 2 -l) -o( /3 cO 3 

(2(8 u; - * £ to 3 ) -1 



. 3 

- °< ft u> 

0 

0 



= 0 



Which, when expanded becomes: 

W(c< ,/3 ,0)) = (3 2 CO 6 - loL ^ 2 (0 4 + 20A) 2 -1)=O 



Thus the equation 

(o/ 2 AJ 6 - 2 «X«) 4 )|G 2 + (+2A) 2 )/S 

gives the constant omega loci except where = 0, 



1 = 0 

0 or 4) = 0. 



(5-32) 



(5-33) 

(5-34) 
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6. Equations of the Constant Bandwidth Curves 



In this section, equations will be developed for displaying curves of 
constant bandwidth upon the parameter plane. 

Let us define what is meant by a constant bandwidth curve: 
n A constant bandwidth curve for G(j £0^) = M is a curve drawn upon the para 
meter plane which specifies the relation between the parameters necessary 
if the transfer function G(s), which is a function of the parameters, is 
to have magnitude M at the real frequency 4) 

Once such curves are obtained for selected values of M, the frequency 
response G(j a) ) may be sketched. Alternately the constant bandwidth curve 
corresponding to some M and specification can be drawn, and from this 

curve values of the parameters can be found which will guarantee that the 
specification is met. 

Consider the following rational transfer function G(s): 



G(s) = 




+ q x 8 + q Q 



Q(s ) 



(6-1) 



11 

r s + r 
n i 



n-1 




n-1 



s 



Where the coefficients of Q(s) and R(s) are of the form: 




(6-2) 



( 6 - 3 ) 



Thus, G(s) may be written: 



m 




( e * + f n£ + g n * + *0 S? 

P P P P 



G(s) = ^ 




( 6 - 4 ) 
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Consider now only real frequency, s =* j , this can be written: 

G( j « ) - 



m 

ZI- 

P-0 


(-l) %P W P (e* +f p (3+g p ^/9 +h p ) + j 


m 

£ 


<-l) 4<p - 1) «J P (e^ 


+£ p e + V^ + V 


even 




odd 






n 

E 

k=0 


(-l)^ k W k (a k e< +b k (S +c k o<^ +d k ) + j 


n 

E 

k=l 


, lx %(k-l) k, 
(-1) 'tv (a k e< 


+ b k e+ c k v / 9 +V 


even 




odd 







(6-5) 



Now let id * 6L?, , some fixed value of 60, and also define: 
b 



n 



- z 

k=0 

even 



(_l)%k . etc. f Qr g , C , D 

b k r r r 



( 6 - 6 ) 



n 



A i = 



2 (-l)*< k_1) 0) k a k ; etc. for B., D. 



E = 
r 



E. « 

l 



k=l 

odd 



m 

z 

P=0 

even 

m 

Z 

p=i 

odd 



(-1) %P 0J P e ; etc. for F ,G , H 
bp’ r r r 



(-l)^ P " 1 ^ a) P etc. for F., G., H. 

b o 111 

p 



(6-7) 



( 6 - 8 ) 



(6-9) 



Using definitions (6-6) thru (6-9), Equation (6-5) can be written: 



(*E r +/3F r + ^G r + H r ) + j(*E. + /3F. + G. + H^) 

G(j A) b ) = (« A r +/ 3 B r + ^3 C r + D r ) + j(*A. +(3B. °i + D i> (6-10) 

# 

Using (6-10) and defining the magnitude of G(j ) as M we obtain: 
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( o< E r + /S F r + o< (S G r +H r ) 2 + ( * E t + /3 F a + o < 0 G^H. ) 2 
M-|G(i6) b )| = (otA+i3B+*/3C+D r ) 2 +(«A i +pB i +<X/3C i +'D i ) 2 

Which when expanded becomes: 

- 0(<*,/S) 

M = Q( * , p ) 

Where: 

0(«, £)=<* 2 p 2 (G 2 + G 2 ) + 2o< 2 /3(E r G r + E.G.) + 2<*/5 2 (F r G r + F.G.) + 

o< 2 (E 2 + E 2 ) +/9 2 (F 2 + F 2 ) + 2^ (E r F r + E.F.) + (6-13) 

2®((E H + E.H ) + 2/6 (F H I- F.ll.) -I- (Il ? ' + H 2 ) 
r r i i 1 r r 1 1 r i 

9( °< > ^ ) = c\ 2 (3 2 (cj c 2 ) •!• 2c^(3 (A r C r !• A.C.) + 2 * (S^B^ + B^.) + 

K*(A J ! A 2 ) I- (J 2 (B 2 + B 2 ) + 2 ^ (A f B r + A.B^ + 

2°<(A r D r + A.D.) + 2|3(B r D r + B.D.) + (D 2 + D 2 ) (6-14) 

Equation (6-12) can be written: 

0 (* >( 6 ) - (M 2 0 (* ,/S)) = 0 (6-15) 

Substituting (6-13) and (6-14) into (6-15) and expressing the resulting 
equation in the form of a quadratic in we obtain: 

(M 2 R 1 -R 2 ) 2 + (M 2 W x - W 2 ) J3 + (M 2 V 1 -V 2 ) = 0 (6-16) 

Where : 

R. = (B 2 + B 2 ) + 2 o< (B C + B.C.) + * 2 (C 2 + C 2 ) 

1 r i' ' r r l ^ v r i 



( 6 - 11 ) 



(6-12) 
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R- = (F 2 + F 2 ) + 2°< (F G + F.G.) + cx, 2 (G 2 + G?) 
2 r i rr 1 1 ' r 1 



W = 2(B D + B.D.) + 2o< (A B +A R +C D+C.D.) + 2c< (AC+A.C.) 
1 rr 11 rriirrii rrii 



W„ = 2(F H + F.H.) + 2c* (E F +E.F.+G H +G.H.) + 2©< (EG +E.G.) 
2 rr ii rriirrii v rrii 

V, =* (D 2 + D 2 ) + 2«((AD + A.D,) + <* 2 (A 2 + A 2 ) 

1 r i rr ii r i 



V„ = (H 2 + H 2 ) + 2<*(E H + E.H.) +o< 2 (E 2 + E 2 ) 

2 r l rr ii r i 



We may now solve the quadratic equation (6-16) for p as a function of , 

U ) and M. This solution is: 

- (m 2 w 1 -W ) + /( -W ) 2 - 4(M 2 R -R )(m 2 v.-v ) 

3 _ 1 V 1 _ _ (6-17) 

1 ? “ 2 

2(M R^Rp 

The constant bandwidth curve corresponding to some desired M = Mq and 
60 ^ = u) bp can be computed from (6-17) by a digital computer. The 
procedure is to fix M and 60 ^ in (6-17) at the desired values, increment o( 
over the range of interest, and at each value of o< to compute the corres- 
ponding values of /3 . The two roots of (6-17) give the two branches of 
the desired constant bandwidth curve. The reader is referred to section 10 
where constant bandwidth curves are computed and discussed for a loaded and 
null adjusted parallel-tee network. 
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7. Equations of the Constant Q Loci for Second Order Systems 

Recent interest in the Q of second order RC networks, accompanied by 
difficulties in applying the conventional definitions of Q, led Morris (4) 
to define Q mathematically as a parameter of the characteristic equation. 
He rearranged the second order characteristic equation into the form (7-1) 
which defines Q: 

7 2 T 

T s + -^-8 +1=0 (7-1) 



Dutta-Roy (9) has pointed out that Q so defined is equivalent to an 
earlier definition proposed by Bolle (10). Bolle's definition for the 
second order characteristic equation: 



is: 



a 2 s + a s + a Q = 0 



Q = 



<a 0 a 2 ) 



(7-2) 

t ^ 

(7-3) 



In this section equations will be obtained which will permit loci of 
constant Q to be displayed upon the parameter plane for variable parameter 
second order characteristic equations of the form (7-2). It is assumed 
that the variable parameters appear in the coefficients of (7-2) in the 
following manner: 

°k = hk* + C k/ 6 + K + f k (7_4) 

k = 0,1, and 2 

Squaring both sides of Equation (7-3) and rearranging we obtain: 



Q 2 a, 



' a 0 a 2 = ° 



(7-5) 



Substituting (7-4) into (7-5) and writing the resulting equation in the 
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form of a quadratic in 



, we obtain: 



KQ 2 d*-d 0 d 2 ),* 2 + (2Q 2 c 1 d l -c 0 d 2 -c 2 d 0 )ol l + (Q 2 c 2 -c 0 c 2 )I (3 2 + 

K^Vl-VVVo^ 2 + < 2<)2b l C l +2Q2d l f r b 2 <: 0- b 0 C 2- d 2 £ 0- d 0 f 2 )o < + 

<2Q2c i E r c 2 f o- c o £ 2 )l P + 

t(Q 2 b 2 -b 0 b 2 ) <=< 2 + ^VrVo-hoV*- « 2£ r f o f 2>' = ° 



(7-6) 



The desired constant Q loci can be obtained from (7-5) by use of a 
digital computer. The procedure is as follows: 

1. Fix Q in (7-5) at a value for which the locus is desired. 

2. Increment over the range of interest. 

3. For each value of 9 solve (7-5) by the quadratic formula. 

4. The two solutions of (7-5), {3 ^ and /3 ^ plotted as ordinates 
on the parameter plane with as abscissa define the two 
branches of the desired constant Q locus. 

5. Repeat steps 1-4 for each value of Q for which a locus of 
constant Q is desired. 

For examples of this procedure the reader is referred to Section 10 
where several loci of constant Q are obtained for a variable parameter, 
loaded, and null-adjusted parallel-tee network. 
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8. Sensitivity Analysis by Parameter Plane Techniques. 



a. Introduction 

The polynomial equation 




k=0 



(8-1) 



where : a fc = f (q^q^ . . . ,q m > 

has n roots. Each of these roots is a function of the parameters q^, 

v it is the objective of this section to show how each of 

these roots vary as the parameters change. 

th 

Let us denote the i complex root pair of (8-1) by: 

r l ' - ?i "i± J Vl\l !- J 2 < 8 - 2 > 

and the j t ^ 1 real root of (8-1) by: 



r. = CT . (8-3) 

J J 

Let us now define two types of sensitivity; macroscopic sensitivity 
and microscopic sensitivity. 

Macroscopic sensitivity is the sensitivity of the roots (8-2) and (8-3) 
to large changes in the parameters q^. These sensitivities may be defined 
by: 

-<*?»)/<* ? r > A <)-=0, Wr 
S i,r = ( & & ^ ? r > A<^=0, i^r (8-4) 

S tr -<*«-,)/( A f r ) ^o ( =0i 0 r 

r v 

S, denotes the change in the damping ratio 4 of the complex root pair 
i ) r * i 

UJ ‘ 

r^ due to a finite change in the parameter q^. ‘ denotes the change in 
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the undamped natural frequency of the complex root pair r\ due to a finite 



change in the parameter q 



Jf* 



denotes the change in the real root 



rj due to a finite change in the parameter q^. 

Microscopic sensitivity is the sensitivity of the roots (8-2) and (8-3) 
to infinitesmal changes in the parameters q^. These may be defined as: 



* 



i.r 



"‘r - iujdxf r 

r 



if * 9 r 



(8-5) 



<fi 



J> r 



The interpretation of these microscopic sensitivities except that the para- 
meter changes are infinitesmal rather than finite. 

b. Obtaining Macroscopic. Sensitivity from the Parameter Plane 

The procedure for obtaining macroscopic sensitivity directly from 
the parameter plane is best demonstrated by an example. 

Example V 

Let us consider some third order polynomial equation: 

a 3 s 3 + a 2 s 2 + a lS + a Q = 0 ; a fc = f fc (o<,^) 



Let us also imagine that we have obtained the parameter plane for this 
polynomial and that it appears as shown in Figure 3; Let us presume that 
we have selected <X= °( q = 3.8 and (3 = /3 = 2.4 in order to place a pair 

of complex roots at = 0.2 and = 0.7 and that this choice results in a 

real root at (T = - 3.0. We wish to determine the macroscopic sensitivities 

X X to u) <r c- 

S . , S. , S , , S , and S. , for an increasd of one unit in 

01 /9 o< (3 o( p 

both ®< and ^9 . From Figure 3 we determine: 
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Fig. 3 Macroscopic Sensitivity from the Parameter Plane. 



c" 

00 






A* 



Ap=o 






^ A<=c 
AW 



A* 

Aoi 






A/®A«t=o 

A <r 



A* 

AO* 



A(i=0 






0.28-0.20 

4.80-3.80 

0 . 12 - 0.20 

3.40-2.40 

0.80-0.70 

4.80-3.80 

0.77-0.70 

3.40-2.40 

2.20- 3.00 
4.80-3.80 

4.20- 3.00 
3.40-2.40 



= 0.08 
= -0.08 

= 0.10 
=-0.07 
= -0.80 
= 1.2 
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More important perhaps, than these numerical values of the macro- 
scopic sensitivities is the qualitative, but enlightening information 



about the sensitivity of the roots which can be obtained by careful in- 
spection of the parameter plane. If the constant zeta, omega, and sigma 
loci have been drawn for constant increments in each of these quantities 
respectively, then one can immediately see that the most sensitive areas 
are those in which the loci are tightly spaced. This is like looking at 
a contour map of terrain- steep areas are highly sensitive areas. 

c. Obtaining Microscopic Sensitivity from the Parameter Plane 
Microscopic sensitivity can be obtained either by the method of 
Kokotovic and Siljak (5) or by an alternate method to be developed here. 
Kokotovic and Siljak derive relations for sensitivity from the basic para- 
meter plane equations. The new method, to be described here, obtains 
sensitivity relations directly from the non-parametric equations of the 
constant zeta and omega loci which were developed in section 5. In those 
cases where the equations of these loci have already been obtained by pre- 
vious calculation, sensitivity relations can be more rapidly developed from 
these loci equations than by the method of Kokotovic and Siljak. As back- 
ground, the method of Kokotovic and Siljak will be discussed first, and 
then the method of obtaining sensitivity relations from the non-parametric 
equations of the constant zeta and omega loci will be introduced. 

(1) Method of Kokotovic and Siljak 
The basic parameter plane equations are: 



n 




k=0 



(8-6) 
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n 



Z (.. 



g = (,-i) k a. ^ ku k ( -T ) = 0 

k=0 C 

Functionally these equations can be written: 

f( 'f , (0 , a k ) - 0 



(8-7) 



( 8 - 8 ) 

(8-9) 



g( \ , A) > & k ) = 0 

Consider that ^ and A) in (8-8) and (8-9) are implicit functions of the 
a^. The a^ are in turn explicit functions of the q^. 

Thus: 

^ f/^ q i = ( 'b J ) q t )+( <3 f/ £%) ( q i )+( ^ f /^ a k )(«) a k /<P q^ 

(8-10) 

d g/ d q. = (d g/dco tO J c) q £ )+( 3 g/ ) ( 3ft / q i )+( g/<? a fc ) (c?a k /c) q £ ) 

( 8 - 11 ) 

But from (8-6) and (8-7): 



t - (1/ 4) t ) X ^’'u. 



k=0 



/Vi<? ■ 



= A. 1 (8-12) 



Li 

c) f/ ^ t (- 



k=0 



dg/d(A ± = (1/ aJj) 



<1 V lc “ 1 .k I / Vy/ \ a , i 

l > *k^ U k-l ( ^ - A 2 



(8-13) 



n 



l <- 



k=0 



* Bj 



5g/ ^' ? n ' B 2 
K— U 



(8-14) 



(8-15) 



Substituting (8-12) thru (8-15) into (8-10) and (8-11) we obtain: 
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n 



^ f / ^ q, = (1/ u) 



,)( duj /d q ) 2 (-] 



..k-l, k TI 

i- T- — „ ' l) ka k"/ U k-l + 

k=0 



n 



k=0 ‘ k ' 1 

n 

2 vX )(_i)k_1 ^X-i =o 

k=0 k r c K L 



(8-16) 



n 



^g/C>q, = (1/ U) .)( d OJ ./<?q r ) 2 (-l) k_1 ka v 0)-V + 



k=0 



k * k 



n 



jd q ) 2 

c k=0 



<-l) k a k *A’ + 



u 

2 (^ a k /^q r )("D k a).\ =0 



k=0 



Let: 



k=0 

n 



:l = (c) a./^q ) (-l) k oj\ 

2 > r k l r c V 



We may now rewrite (8-16) and (8-17) as: 



(8-17) 



u 

- Z (^/^r) (-l) k "X\-l (8_19) 



(8-20) 



- C l.r 


(8-21) 


<r 


(8-22) 



Equations (8-21) and (8-22) are simultaneous equations, from which the 
desired sensitivities ( d id q^) and (d (O jd q r ) can be obtained. 
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d. Sensitivity Relations from the Constant Zeta and Omega Loci 
In Section 5 it has been shown how the equations of the constant zeta 
loci z ( a k> ) = 0 and the constant omega loci W(a k> u) ^) = 0 can be 
obtained in non-parametric form. From these equations we can obtain the 
relations for the microscopic sensitivity directly: 

Z( V £■ ) - 0 (8-23) 

W(a k> to- ) = 0 (8-24) 

from which: 

i/c5» a k = -(^Z/^a k )/(c) Z/ .) (8-25) 

providing ^ Z/ d) £ ^ ^ 0 “ 

J d* = -(<P W /<J a )/(^W/ .) (8-26) 

1. K K 1 

providing 3 id ^ £ 0 

For an example of sensitivity computed from the non-parametric equa- 
tions of the constant zeta and omega loci the reader is referred to 
Section 11 where root sensitivity to parameter variations in a lattice net- 
work is discussed in detail. 
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9. Parameter Plane Analysis and Design of the Symmetrical Parallel- 
Tee Network 

a. Introduction 

The symmetrical parallel-tee network shown in Figure 4 has 
been the subject of several papers (9,10,11). The reason for this interest 
stems primarily from the fact that it is possible to obtain complex zeros 
of transmission using the network. These may be placed anywhere in the 
left half s-plane and as far as 30 degrees from the imaginary axis into 
the right half s-plane by suitably choosing the network elements. From 
the transfer function of the network. Equation (9-1), it can be seen that 
the poles and zeros of transmission are functions of the parameters o( and 

0 • 



E o< p) = <*/3p 3 + 2/Qp 2 +2/3p + l (9-1) 

E i (p) Cffi P 3 + (tf+2/S + 2* ) P 2 +(*+2f3+2)p+ 1 

where: p = RCs 

Barker and Rosenstein (12), using root locus techniques, investigated 
this network. Their analysis, while complete and well done, was compli- 
cated considerably by the inherent single parameter capability of the 
root locus technique. This section approaches the same problem with a 
two parameter method of analysis and design - the parameter plane. 

By employing parameter plane techniques, graphs are prepared for 
the polynomials which form the numerator and denominator of the transfer 
function. These graphs, with as abscissa and ^3 as ordinate, display 
directly the roots of the polynomial which result for any choice of the 
parameters o< and ^3 . From these parameter planes for the network, the 
regions of the s-plane in which poles and zeros of transmission are possible 
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are immediately evident. Likewise values of the parameters o< and ft 
necessary to place either poles or zeros at some desired point within these 
regions can be directly read from the curves. The parameter plane is there- 
fore a tool of both analysis and design. 

Loci of constant bandwidth are drawn upon a parameter plane for this 
network. These consist of a family of curves along any one of which the 
frequency at which the response of the network is down N db is constant. 

For any choice of o< and (3 the frequencies at which the response is down 
N db can be read directly from these loci. Loci for the conventional -3db 
bandwidth are presented in this section. 

Thus it is possible to select values of the parameters o< and 
which guarantee suitable s-plane locations for certain poles and zeros 
while indicating the bandwidth determined by this choice. Conversely o< 
and p may be chosen to guarantee upper and lower -3db frequencies while 
offering a selection of poles and zeros which may be read from the curves. 

The curves which are presented in this section are frequency normalized 
and hence universal. They may be used to investigate any symmetrical paral- 
lel-tee network of the form shown in Figure 4. 

The advantage of the parameter plane over other methods which have 
been used to analyze and design this network are the rapidity and ease with 
which the capabilities and characteristics can be determined from the curves 
and with which the parameters may be selected in order to produce desired 
network behavior. 

b. Approach to the Problem 

The objective of the analysis to follow is to determine 
in what manner the poles, zeros, and -3db frequencies of (9-1) change as the 
parameters o( and ^3 are varied. This can be accomplished by preparing three 
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parameter planes; one each for the poles, zeros and -3db bandwidth of (9-1). 
Since some of the equations of these parameter planes have already been 
developed in earlier sections, they will be recalled as appropriate. 



c . Poles of the Transfer Function 

Since the network shown in Figure 4 is passive RC, all of 
the poles will lie along the negative real axis in the p-plane. Let <T 
denote the value of p at which such a pole of (9-1) occurs. We may thus 
assert that poles occur at those values of <T which satisfy (9-2): 

* £<r 3 + (2 p, + 2o(@ +o()<r 2 + (2 (5 + 2* + 2)0*+ 1-0 ( 9 - 2 ) 

Rewriting (9-2) in the form of a plane algebraic curve in the 0 

plane we obtain: 

< <r 3 + 2 <r 2 ) o (|3 + (2 <r 2 + 2 (T > (i + (<r 2 +cr )c< + (2<r+ i> = o 

(9-3) 

Equation (9-3) represents a family of hyperbolas in the o< - plane. One 
hyperbola occurs for each selected value of (f . Since (9-2) is satis- 
fied by three values of 0* for each choice of and f$ , it follows that 
each point in the parameter plane will have three hyperbolic branches pass- 
ing through it. 

The curves represented by (9-3) can be hand calculated for selected 
values of (T* , but since this is a tedious process, a digital computer was 
used to obtain them. The resulting curves which specify the pole values 
for any choice of the parameters are shown in Figure 5. One enters this 
figure with values of o< and and then reads the values of the three 
poles directly. Alternately one may select values of o( and ft to produce 
some selected pole values. 
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d. 



Zeros of the Transfer Function 



In this instance the problem is much more interesting since 
complex zeros are possible. The zeros are the roots of: 

cj p p 3 + 2/Sp 2 +2/S P+1 = 0 (9-4) 

Since o< and fi are presumed real it can be asserted that (9-4) has at 
least one real root. As before let p = (X in (9-4) and rewrite the re- 
sulting equation as a plane algebraic curve in the oc - |3 plane: 

0-3 + (2 <r 2 + 2<r ) ^3 +1 = 0 (9-5) 

Equation (9-5) also represents a family of hyperbolas in the o< - / 3 plane. 
These are shown in Figure 6. Note that in some regions it is possible to 
have three real zeros. 

In describing the location of complex zeros of (9-4) the zeta and omega 
notation familiar to control systems theory will be used: 

s = - ^ to + j CO \j 1- ^ 2 ' (9-6) 

The non-parametric equations for the constant zeta and omega loci which were 
developed for this network in Section 5 will be employed. The equations of 
these loci are repeated below. For their derivation the reader is referred 
to Section 5. 

z( * , /S , 7} ) = - */3 3 (4 j 2 -l) 3 [(-32 * f 2 + 16)/3 2 + (64 * f 4 + 

16 o( 7j 2 -32 £ 2 -8*)^6 + (-64* 2 J 6 + 48* 2 JT 4 - 12* 2 ^ 2 + •< 2 ) =0 

(9-7) 

w , /3 , 00 ) = °< OJ 3 (* 2 /S 2 (J - 2 c* £ 2 4) 4 + 2 p od 2 - 1) = 0 

(9-8) 
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A digital computer, programmed as described in Section 5, was used to 
obtain the constant zeta and omega loci which are shown in Figure 6. 



e . -3db Constant Bandwidth Curves 

One may obtain the desired -3db constant bandwidth curve by 
employing the procedures which were developed in Section 6. Referring to 
that section we note that several quantities must be substituted into 
Equation (6-17) in order to obtain the curve which we desire. In order to 
facilitate substitution into (6-17) the following quantities are tabulated: 
From the transfer function (9-1) using the notation of Equations (6-2) 



and (6-3): 



k 


fk 


fk 


fk 


fk 


fk 


fk 


fk 


fk 


0 


0 


0 


0 


1 


0 


0 


0 


1 


1 


jW b 


j2^ 


0 


j2*) b 


0 


J^b 


0 


0 


2 




-2*0 2 
b 


b 


0 


0 


- 2*1 

b 


0 


0 


3 


0 


0 


-&1 


0 


0 


0 


- 5*1 


0 


From Equations 


; (6-6) 


thru 


(6-9) - 


using values 


from (9-9) 


one 


obtains: 



A 

2 

Real -60. 

b 

Imag £) b 

From Equation 

R = 
1 



B 


C 


D 


E 


F 


G 






1 


0 


*V 

CM 

1 


0 


24) 

b 


< 


2 *b 


0 


2* 

b 


-4, 3 

b 


(6-16) 


using 


the values 


from 


(9-10) one 


obtains 


4 fc) 4 
b 


+ 

b 


+ 4o< <J 4 + 
b 


4« 2 


a! + « 2 of 

b b 





4 2 2 6 4 

4AJ b + 4i °b + ^ < -4*^b 

2 * 2 a ; 4 - 4 a)* 

b b 



H 

1 

0 



(9-10) 



(9-11) 
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W„ = -4 60 2 
2 b 

V = 1 + °< 2 6t)^ + <x 2 tol + 2<Xtd 2 + l\0d 2 
1 b b b b 




The quantities (9-11) are used to digitally compute values of Equation (6-17) 
corresponding to selected values of 40 ^ with M = 1/ \J 2^as required for the 
-3db curves. The resulting curves are displayed in Figure 7. 

f . Interpretation of the Parameter Planes 

The parameter plane for the poles of the transfer function 
which is shown in Figure 5 reveals that for 0 < o( , fa 2.5 the three 
poles are confined to regions along the negative real axis as shown below: 




The parameter plane shown in Figure 6 for the zeros of the transfer 
function reveals the following information regarding the zeros: 

(1) The positive axis coincides with the $ = -% locus 
indicating that zeros may be placed as far as 30 degrees from the imaginary 
axis into the right half p-plane. 

(2) The zeta = 1 locus separates the region of 3 real zeros 
from the region of 2 complex conjugate zeros and 1 real zero. 

(3) The zeta = 0 locus is a straight line with equation 

= 4 ft . This indicates that if one is interested only in those cases 
where a pair of conjugate zeros are placed on the j 40 axis, then he may 
substitute this linear relation between the parameters into the network 
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^uations thus reducing the problem of analysis and synthesis to one of a 
single parameter. 

(4) The zeta = \ locus coincides with the omega » 1 locus 
indicating that an additional degree of freedom exists which may be used 
to specify other poles or zeros in addition to the specification that a 
pair of complex conjugate zeros occur at f = %, ^ = i. 
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Fig. 4 Symmetrical Parallel-Tee Network 
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CURVE LABELS REFER TO 
NEGATIVE REAL ROOT VALUES 

Fig. 5 Poles of the Symmetrical Parallel-Tee Network 



50 



A 




a 

Fig. 6 Zeros of the Symmetrical Parallel-Tee Network 
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Fig. 7 -3db Constant Bandwidth Curves for the Symmetrical Parallel- 

Tee Network 
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10. Loaded and Null Adjusted Symmetrical Parallel-Tee Network - 

Parameter Plane Analysis 
a. Introduction 

The null - adjusted symmetrical parallel-tee network provides 
a means of obtaining frequency selectivity without the use of inductors. 
Hence it is a very useful network at low frequency where inductor size be- 
comes unacceptably large as well as in integrated circuits where the current 
technology does not permit inductors of reasonable size to be integrated. 

A number of articles have considered the zero source - infinite load 
resistance case (9,13,14), and a few have considered the effects of other 
finite combinations (10,11,12). These latter articles generally consider 
either the source or resistance to be fixed and the load resistance to be 
variable or vice-versa: 

This section considers the problem from a different point of view and 
with a two parameter tool of analysis and design - the parameter plane. 

•V 

The two parameter nature of the problem is evident by inspection of Figure 
8 and Equation (10-1): 



E „<p> = 



£ <p 2 ± i) 



(10-1) 



E l (p> (4*|« +o( +/S)p 2 +(4 ®f/9 +4 of +4(9+2)p + (=< +/S+ 2) 



where: p = RCs, oc = R^/R > = R^/R 

Parameter planes will be obtained for this transfer function. They 
will show the poles, Q, and frequency response characteristics associated 
with this network which result for any choice of the parameters o( and ^9 with- 
in the range 0 <, , /3 < 3. Alternately and /S may be chosen from the 

curves to obtain desired characteristics within the capabilities of the 
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network. 



Loci of constant Q discussed in Section 7 are also obtained. These 
loci show directly what value of Q results for any selection of parameters 
c< and (3 . Conditions are specified which permit Q to be related to the 
bandwidth of the network. 

Constant bandwidth curves are obtained as described in Section 6. 

For any choice of parameters these curves immediately give the upper and 
lower frequencies where the response of the network is down N db from the 
response at infinite frequency. Curves are presented for several values 
of N, thus permitting the frequency response to be sketched for any chosen 
values of source and load resistance. 

b . Poles and Zeros of the Transfer Function 

The parameter plane graphically displays the roots of a 
polynomial whose coefficients are functions of two variable parameters. 

In the general case these polynomials (which may be the numerator or denom- 
inator of some network function) can have complex as well as real roots. 

In the case of Equation (10-1), however, the following observations are 
made : 

2 

The zeros are fixed at p = - 1 

The poles are all real and lie on the negative real axis 
in the p-plane since the network is passive RC. 

Thus it is necessary to investigate only the poles of (10-1) and this is 
considerably simplified by the fact that the poles are all real. To obtain 
the parameter plane for these poles let p = C T be a pole of (10-1). Making 
this substitution into the denominator of (10-1) and rearranging the result 
ing equation into the form of a conic section in the ot - /$ plane one obtain 
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(4 cr 2 + 4cr ) «*^3 + ( cr 2 + 4<r + l)<* + (<r 2 + 4cr +i)/3 + (2C/+ 2) = 0 

( 10 - 2 ) 

This equation represents a family of hyperbolas; one hyperbola for each 
value of CT . These are plotted for several values of CT and are shown in 
Figure 9. One may read the poles, which result for any choice of ^ and, 

^ , directly from the curves. 



Example VI 

Assume that a parallel-tee notch network has been designed to reject 
60 cps. The notch filter has been designed so that R = 2660 ohms and C = 1 
microfarad. The network is driven from a source resistance of 5000 ohms 
and drives a load resistance of 2500 ohms. It is required that the poles 
of the transfer function be obtained. From Figure 8 one obtains: 



o< = 5000 = 1.88 
2660 

0 = 2500 =0.94 

2660 



From Figure 9 the normalized poles are read: 

cr; = -1.78 
O'-L = -0.27 



Which, when un-normalized become: 



= 



-1.78 



2.66 x 10 
-0.27 



-3 



u 2 2.66 x 10 

The network transfer function is thus: 



-3 



= -670 
= -102 
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E o< s) = /3 f s 2 + Jo 2 ] 

E i (s) 4 ctfi +o< + /3 (s +(T{) (s +<T' 2 ) 

= 0.95 f s 2 + 377 2 -I 

(s+670) (s+102) J 

The curves of Figure 9 can be used to place one of the poles at a de- 
sired location while maintaining some control over the location of the other 
pole. This is accomplished by changing , or R. Changing R and 

assuming that C is also changed appropriately to keep the zeros fixed changes 
both o< and /S . 

Example VII 

Suppose that in the previous example one is able to change the source 
resistance and desires to do so in such a manner as to cause the pole at 
s = -670 to move to s = -750. Thus the new normalized pole location be- 
comes: p^ = -2.00. Since the load resistance is fixed at 2500 ohms, /S> 
remains 0.94. From Figure 9 note that the CT = -2.00 locus intersects the 
P =0.94 coordinate at o< = 1.07, and that the other pole is thus located 
at p^ - -0.335. When un-normalized these pole locations become: 

s x = -750 
s 2 = “126 

The new transfer function becomes: 

E o< s) = 0.156 [ S 2 + 377 2 ] 

E i (s) (s + 750) (s + 126) 

The new source resistance required becomes: 

R g = ©< R = 1.07(2660) = 2840 ohms 



56 



Constant Bandwidth Curves 



c . 

By making the substitution p = j od as discussed in Section 
6, the real frequency transfer function is obtained for the loaded and null 
adjusted symmetrical parallel-tee network: 



Where: 



E o<*> > 

E t ( CO ) 



K ( « , jS ) (G(«,/8,6))] 



6 



(10-3) 



(10-4) 



G(«,/S ,*)) = (4^/6 +°( ±j3 )(l-4) 2 ) (10-5) 

’ [(* + /S+2)-(4*/3+o(+/S ) CO 2 ] +j[ (4 +4 +4^6 +2)] 

By inspection G(®( , ^9 , CO )— >1.0 as (O -♦ °o . At any selected frequency 
<0 , and attenuation level M at that frequency, the relation 

G( of , 0 b ) = M (10-6) 



defines a constant bandwidth curve in the - /3 plane. This curve specifies 
the source-load resistance relation requires for (10-6) to be satisfied. 

The attenuation level M can be thought of as the attenuation from the in- 
finite frequency value of E q (00 )/E^(co )> which is K( , {S ): 



K( o( , ) = V °° > 

E i (00 ) (10-7) 

In order to obtain the equation of the constant bandwidth curves Equation 
(10-5) is substituted into Equation (10-6) and both sides of the resulting 
equation is squared. This squared equation can be solved for /3 as a 
function of o( , CO and M: 
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(10-8) 






-(M^-Wp + f(M 2 W 1 -W 2 ) 2 -4(M 2 R 1 -R 2 )(M 2 V 1 -V 2 )]' 
2(M 2 R 1 -R 2 ) 



Where in the notation of Section 6: 

R 1 = <0 4 + 14A) 2 + 1 + 4o< (2 CO 4 + 6 A) 2 ) + 16<* 2 ( A) 4 + A) 2 ) 

r 2 = A) 4 - 2 A) 2 + 1 + 4<* (2A> 4 -4 4) 2 +2) +16* 2 ( aJ 4 -2 U} 2 + 1) 

W x = 12A) 2 + 4 + 2<tf (AJ 4 + 14 A) 2 + 1) + 4* 2 (2A) 4 + 6*) 2 ) 

W 2 - 2 *< ( A) 4 -2 AJ 2 + 1) + 8 °< 2 ( AJ 4 -2 A) 2 + 1) 

V]L = 4 + 4 A) 2 + 4*(3A) 2 + 1) +c* 2 (a) 4 + 14A) 2 + 1) 

2 4 2 

V 2 * o< ( (o -20J + 1) 

Assisted by a digital computer, Equation (10-8) was used to obtain the 
constant bandwidth curves for values of M corresponding to -3db, -6db, and 
-20db. These curves are shown in Figures 10, 11, and 12 respectively. 
These curves may be used to sketch the frequency response as indicated in 
the following example. 



Example VIII 

For the conditions specified in Example VI, namely = 1.88 and P 
0.94, the curves of Figures 10, 11, 12 are used to obtain: 



Db down from infin- 
ite frequency gain 


s 


cd h 
1 2 


-3 


0.50 


2.65 


-6 


0.61 


1.84 


-20 


0.901 


1.115 
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The infinite frequency gain is: 

K(1 .88, 0.94) = 0.095 
The zero frequency (DC) is: 

E o<°> = @ = 0.195 

E i (0) cX+^H- 2 

The attenuation at O) = 1.0 is infinite since the network of Example VI 
was designed for this condition. 

The above information may be used to sketch the frequency response 
of the network. 

d. Constant Q Loci 

As pointed out in Section 7, the Q of a second order network 
may be defined as: 

- -th , n nd . 

^ _ geometric mean of 0 and 2 degree terms 

l St degree term 

where these terms refer to coefficients in the denominator or the network 
transfer function. For the loaded and null-adjusted symmetrical parallel- 
tee network, Q becomes: 

O = +*+/ 3 )(* + /3 + 2)]^ 

^ 4 + 4* +4/3+2 (10-9) 

Squaring both sides of Equation (10-9) and rearranging the resulting equa- 
tion into the form of a quadratic in (3 obtains: 

[16Q 2 * 2 +(32Q 2 -4)o< +(16Q 2 -1) ] £ 2 + 

[ (32Q 2 -4) ®( 2 + (48Q 2 -10)o( + ( 16 q 2 - 2 )]/S + 

[(16Q 2 -1)o< 2 + (16Q 2 -2)o< + 4Q 2 ] =0 (10-10) 
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Equation (10-10) was used to obtain the constant Q loci displayed in Figure 
13. 

Q defined in this manner is a mathematical parameter of the network as 
discussed in Section 7. It is not in general directly related to the 
frequency response, however, if the infinite frequency response is made 
equal to the DC response as explained later, then the following relationship 
holds: 

Q = f ° (10-11) 

2 “ 1 

Where: f^ is the notch frequency 

f^ is the lower -3db frequency 
f 2 is the upper -3db frequency 

It has been shown (10) that the null adjusted symmetrical parallel- 
tee network under the conditions R g = 0, = has a Q equal to 0.25* 

Figure 13 shows that Q may be made larger than 0.25 by adding source re- 
sistance. The maximum attainable Q is .353* and occurs at R g = R^ = 0.707R. 
e. Locus of Symmetrical Frequency Response 

Also displayed on Figure 13 is the locus of symmetrical 
frequency response. This locus defines the source - load condition re- 
quired for the DC response to equal the infinite frequency response. 

2 

Cowles (10) has shown that this requires R R - R /£ or equivalently 

S L 

(6 = %. Figure 13 shows that this locus also passes thru the point 
of maximum Q. Thus a loaded and null adjusted symmetrical parallel-tee 
network whose source and load resistances are chosen so as to maximize 
Q will also have a symmetrical frequency response. As pointed out above 
in sub-section d, the network under these conditions also is the most 
selective by virtue of its optimum Q. 
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Fig. 8 Loaded and Null Adjusted Symmetrical Parallel-Tee Network 
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Fig. 9 Poles of the Loaded and Null Adjusted Symmetrical Parallel 

Tee Network 
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Fig. 10 -3db Constant Bandwidth Curves for the Loaded and Null 
Adjusted Symmetrical Parallel-Tee Network 
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Fig. 11 -6db Constant Bandwidth Curves for the Loaded and Null , 
Adjusted Symmetrical Parallel-Tee Network 



/ 
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Fig. 12 -20 db Constant Bandwidth Cui'ves for the Loaded and Null 

Adjusted Symmetrical Parallel-Tee Network 
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CONSTANT Q LOCI 




Fig. 13 Constant Q Loci and Locus of Symmetrical Frequency Response 
for the Loaded and Null Adjusted Symmetrical Parallel-Tee 
Network 
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11. RC Lattice Network - Complex Zero Synthesis and Sensitivity Analysis 
by Parameter Plane Techniques 
a. Introduction 

The RC lattice network which is investigated is shown in 
Figure 14. The parameters o< and f$ are real and variable in the arbitrarily 
selected range 0 < 2.5. The transfer function of this network 

is given by: 



e d (p) 

E.(p) 



oQQ p 2 + ( - <S + pp + i (li-D 

<*/3 p + ( o( (3 + (3 + i)p + i 



where : p = RCs 

Frequency normalization makes the results which are obtained universally 
applicable . 

Parameter planes which are obtained graphically display the poles and 
zeros of (11-1) which occur for any choice of parameters and /5 within 
the range specified. These parameter planes show in what manner the poles 
and zeros change as the parameters are varied. Using the parameter planes 
one may select parameter values which are required in order to place poles 
or zeros at desired p-plane locations within the capabilities of the net- 
work, The analysis performed shows that it is possible to readily place 
zeros anywhere in the complex region of the p-plane. 

Macroscopic root sensitivity can be determined by inspection of the 
parameter planes. Equations are obtained as described in Section 8 which 
permit microscopic root sensitivity to be computed for selected root loca- 
tions . 
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b. 



Poles and Zeros of the Transfer Function 



The poles of the transfer function are given by the roots 



of : 



<*y 9 p 2 + ( *(/3 + /3 + i)p + 1 « o 



(11-2) 



Since the network of Figure 14 is passive RC , the roots of (11-2) all lie 
on the negative real axis in the p-plane. Let p = (T be such a root of 
(11-2). Making this substitution into (11-2) and rearranging the result- 
ing equation into the form of a conic section in the o< -^9 plane one obtains: 



(Q * 2+<r >° < / 3 + <T ( 3 + (<r + i) = o (n-3) 

Equation (11-3) represents a family of hyperbolas in the o( - plane; 
one hyperbola results for each value of (T . This equation is solved for 
as a function of o< and <T and used to' graph the curves displayed in 
Figure 15: 

* = + i) (n-4) 

( <r + <r ) * + <r 

The zeros of the transfer function are the roots of: 



p 2 + 




+ 1- /9 )p + i = o 



(n-5) 



Following the procedures of Section 5 the constant zeta and omega loci are 
obtained in non-parametric form: 

(1) Constant Zeta Loci 



From which: 



f( to ) - 2 f # (3 V 2 + - ofrs -1)4) + o 

g(tO ) = 0U) 2 + ((S - -1)£0 + 2 f 



0 

0 






2 #3 (p-'.'c-l) 



7 i 
0 



2p K- 



0 (£ -72-1) 



0 0 

6-72-1) 0 

2? 0 

(3-76-1) 2 X 



( 11 - 6 ) 

(11-7) 



( 11 - 8 ) 
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Expanding (11-8) one obtains: 
.2 . 2 



Z(* ,{S [((X -2c<+l)(S z + (-4 +2* -2)|fl + 1] = 0 (11-9) 



Thus: 



(4f 2 o/ -2<X+2) ± / (-4 y 2 «( +2^ -2) 2 - Aft* 2 -2 <k'+1 ) 



( 11 - 10 ) 



2(o< -2o( +1) 

2 

gives the constant zeta loci everywhere except where 4 J = 0, 

Equation (11-10) was used to obtain the constant zeta loci shown in Figure 
16. 

(2) Constant Omega Loci 

Using Theorem II, Appendix II, for the cases j = 0 and j = 1, one 
obtains: 



-( fi - «/3 -l)«J + u) 2 




1 

Cj 

»-» 

1 


(<Xp 2 -l ) 0 




1 

CM 



= 0 



(11-11) 



From which: 



( 11 - 12 ) 



w( ex' , ,4J ) = <*ys do ( o('^3 od -1) = 0 

2 

Thus o f{& = 1/ OJ gives the constant omega loci everywhere except where 

2 

of ft (aO - 0. Equation (11-12) was used to obtain the constant omega loci 
shown in Figure 16, 

c . Root Sensitivity 

(1) Macroscopic Sensitivity 

Root variations due to large parameter variations are read- 
ily observed by inspecting Figures 15 and 16. In the case of the sensitivity 
of the zeros of the transfer function, one can see that zeros with 0.3 < 

< 0.7, and 0.5 < CO K 1.0 are relatively insensitive to parameter varia- 
tions. Given a specification on permissible zeta and omega variation, one 
can determine, practically by inspection, the required of and tolerances 
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which will meet the specification. 

(2) Microscopic Sensitivity 

Zero sensitivity due to infinitesmal variations in the para- 
meters may be readily computed from formulas derived from the non-parametric 
equations of the constant zeta and omega loci as discussed in Section 8. 
Substituting derivatives obtained from (11-9) and (11-12) into (8-25) and 
(8-26) one obtains: 

c>^ i = - u) (11-13) 

d* 2 * 

dlO i = - 60 (11-14) 

dp 2 $ 

where: 2 +2)N 

i = - ? (8 ? 2 * & -3 2 Q 2 + 4 * (S 2 -4 * (3 - /3 2 +2 /3 -1 ) (11-15) 

2 c< (8 $ 2 *yS - <* 2 ($ 2 + 2 ‘K'/S 2 - /Q 2 -2 +2/6 -1) 

j £ i a - y (8 ? 2 o<VS -3 * 2 2 + 6 2 -3 2 -4 */3 +4 /6 -1 ) (11-16) 

2/3(8/ V/tf - * 2 0 2 +2 w'/S 2 - /3 2 -,2 V/3 +2 /S -1) 

where: /3 (8/ 2 c*/S - * 2 /S 2 + 2 V/S 2 - (3 2 -2 V/3 +2/S -1) 4 0 

d. Interpretation of the Curves 

Inspection of the parameter plane for the zeros of the 

transfer function shows that it is possible to place a pair of complex con- 

jugate zeros anywhere in the s-plane by appropriately choosing the para- 
meter values &( , & and the frequency normalizing factor RC. 
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Fig. 14 RC Lattice Network with Variable Parameters 
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Fig. If Poles of the Variable Parameter Lattice Network 
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Fig.j.6 Zeros of the Variable Parameter Lattice Network 
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12. Application of the Method to other Problems 

This section contains examples of several networks whose transfer 
functions are examined to determine the applicability of parameter plane 
methods . 

a . Pi Network 

The pi network shown in Figure 17 is useful for matching a 
resistive source to a resistive load at a fixed frequency. The source and 
load resistances are chosen as variable. 










Fig. 17 Pi Network with Variable Source and Load Resistances 
The transfer function of this network is given by: 

V s > = (S 

E i/ S ^ °*yS CjC^s 3 + (a<'C 1 L +(9C 2 L)s 2 + ( °< (S C + or '(3 L)s +(o<+(6) 

( 12 - 1 ) 

Equation (12-1) shows that the parameters appear in the coefficients in a 
manner that permits a parameter plane to be constructed which will show 
how the network poles are determined by the values of the variable para- 
meters . 

b . Unsymmetrical Parallel-T Network with Variable Parameters 

The network shown in Figure 18 is of particular interest here 
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because of the way that the parameters appear in the coefficients of the 
transfer function. 
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Fig. 18 Unsyrametrical Parallel-T Network with Variable Parameters 
The transfer function of this network is given by: 

E o (s) = o /|9 P 3 +( <* /3 2 + < 3 2 )p 2 -K 2 + /9 2 )p + */3 ( 12 - 2 ) 

E i (s) of/3 p 3 +( of/S 2 +/5 2 +o<(S + °<+j5)p 2 +( <X/3 2 -i- /3 2 + ofp+^ + ^p + 1 

2 

First inspection of (12-2) indicates that the (6 terms will complicate the 
application of parameter plane techniques. The following is noted however: 

(1) The network is passive RC, therefore the poles are all real and lie on 
the negative real axis in the p-plane. Hence p = (T* may be substituted into 
the denominator of (12-2). This denominator is then equated to zero and 
the resulting equation solved by the quadratic formula for as a function 
of o( and . Thus the parameter plane for the poles may be constructed. 

(2) The quantity (S may be factored from the numerator of (12-2). The 
resulting polynomial is of the proper form to be analyzed by parameter plane 
methods . 
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c . Bridged-T Network with Variable Parameters 

The bridged-t network shown in Figure 19 can be used to gener- 
ate complex conjugate zeros in the left half s-plane. 

fin 




Fig. 19 Bridged-T Network with Variable Parameters 
The transfer function of this network is: 



E o (p > = o/ /9 s 2 4- (o<+ l)s + 1 

2 



E.(p) 



(12-3) 



<*(3 s + (o< +/?+ 1 ) s + 1 

where p = RCs 

The form of Equation (12-3) readily permits parameter plane analysis and 
design . 

d. Bridged Twin-T Network with Variable Parameters 

The network shown in Figure 20 is an RC Bridged Twin-T 
Network with variable bridging and shunting elements. 




S 

Fig. 20 Bridged Twin-T Network with Variable Parameters 
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The transfer function of this network is given by (12-4): 



E^P) 

Zo</8<f+2/8rM<t)P%( i fiJ+2wJ-&* t + 2-< f )p 2 +(2fi - ^£+2&)Pi-l 
(o<fi)f* < f)p i 4-(2d/3t t J+2/2>x'cf £+*( sar* +yj) p S i- 

2 

(x/S 4-2J + 2) p +/ 



Where p ® RCs 

Inspection of (12-4) reveals that one may construct several parameter planes 
depending on which pair of the quantities o/ , > S > are considered 

variable and which pair are considered fixed. The following possibilities 
are available: 



Variable Parameters 



Fixed Parameters 



e . 



> ft 


t,6 


t ,<r 


<* , /3 


o< 




/s ,s 


o< , r* 


o( > £ 


/6.V 


j3>Y 


» <f 



Capacitive Divider Matching Network 



The network shown in Figure 21 is frequently used for 



frequency selective coupling between high output resistance and low input 



resistance amplifier stages. 

o< 
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The transfer function of this network is: 



E o (s) = /9y LS 2 (12-5) 

E i^ ( l 2 )s 4 +( /S>^J'l 2 +^l 2 )s 3 +(o(^9^ l-lV + ^J’ L4-/9Vl)s 2 + 

(oty 1 L+*/3<T + <v/3V + l) s + <* 

Equation (12-5) shows that it is possible to construct parameter planes 
with the capacitances^ and J" variable and the source and load resistances 
o( and fixed or vice versa. 

f . Two Section Resistively Loaded LC Low Pass Filter 

The LC low pass filter shown in Figure 22 has variable 
source and load resistances. 




Fig. 22 Two Section Resistively Loaded LC Low Pass Filter 
The transfer function of this network is: 



E.(s) = 

Viy 



0 



(/S L 1 L 2 C 1 C 2 )s 4 + ( o(/S c 1 c 2 l 2 +l i l 2 c 1 )s 3 + 



(12-6) 



(IBl c^IBl^^ £ l 2 c 2 )s 2 + ( C^^C^+L^s + 

(«V+/8 ) 

The form of this function permits application of parameter plane techniques. 
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g. Butterworth Filter Pole Sensitivity to Coefficient Variations 
Parameter plane techniques can be used to investigate the 
effects upon pole locations of Butterworth Filters due to variations in the 
coefficients of the transfer function caused by component tolerances, aging, 
etc. As an example, consider the transfer function of a normalized 4th 
order Butterworth filter: 



E q (p) 



1 



(12-7) 



E i (p) p 4 + 2. 613p 3 + 3.414p 2 + 2.613p + 1 

Assume that the last two coefficients in the denominator of (12-7) are 
variable. Denote these variable coefficients as o( and fS and equate the 
denominator of (12-7) to zero. This results in the polynomial equation: 



4 

P 



+ 2.613p^ + 3.414p 2 + p + 



= 0 



( 12 - 8 ) 



The parameter plane for this polynomial is displayed in Figure 23. When 
°( = 2.613 and ^3 = 1.0 as is the normal case, the poles of (12-7) are 
equally spaced on the unit semicircle in the left half p-plane. By enter- 
ing Figure 23 with values of & and /3 , the resulting poles of (12-7) 
may be read directly. 
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13. Algebraic Design of Electric Networks by Parameter Plane Techniques 



a. Introduction 

Parameter values which will guarantee that certain poles and zeros 
of a rational network function appear at desired s-plane locations and at 
the same time guarantee that other specifications are met may be algebraical- 
ly determined by employing parameter plane techniques. This may be accomplish- 
ed without the necessity of drawing parameter plane curves and thus, more 
than two adjustable parameters may be simultaneously considered. In general 
it is necessary to have at least n adjustable parameters in order to make n 
specifications. While the design technique does not mathematically demand 
it, some preliminary investigation is desirable to determine whether or not 
a physically realizable set of parameter values is likely to be found which 
will meet the specifications. If such an investigation is not performed 
then it may be the case that the design technique will generate only sets 
of parameter values which are unrealizable. 

b. Types of Specifications Which May be Used 

Any specification on the network function or on the network be- 
havior which can be expressed as an algebraic equation in the adjustable 
parameters can be used in conjunction with this technique of design. These 
specifications can include both pole-zero locations and performance specifi- 
cations. Some of these are: 

(a) Pole-zero locations. The parameter plane equations of Section 
4 can be used to specify desired pole-zero locations. 

(b) Zeta of a pair of complex roots may be specified without 
specifying the corresponding omega, or vice versa, by using the non-para- 
metric equations of the constant zeta and omega loci presented in Section 
5. 
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(c) Bandwidth. The equations of the constant bandwidth curves 



developed in Section 6 can be used to specify bandwidth. 

(d) Q. The equations of the constant Q loci, Section 7, may 
be used to specify Q for second order systems. 

(e) Maxima and minima in the frequency response can be specified 
by equating the derivative of the magnitude function to zero at the fre- 
quency where the maximum or minimum is desired. 

(f) Infinite frequency or DC gain of a network may be specified 
by expressing these quantities as algebraic equations in the adjustable 
parameters . 

c. The Technique of Algebraic Design 

(1) Pole Specifications 

Consider either the numerator or denominator of a rational 



network function written as the polynomial equation: 



n 

k=0 



a s^ = 0 
k 



03-1) 



where the coefficients a^ are functions of the adjustable parameters q^: 



= * ^ 2 * * * * 9 ^i^ (13-2) 

Parameter plane theory provides a number of algebraic equations which 
relate the root factors CO 9 <J^ to the parameter values q^. These are: 

(a) The basic parameter plane equations for real and 



complex roots: 



1. For real roots: 



n 

z 

v=n 



a k <r 



k 



0 



(13-3) 



2. For complex roots: 




= 0 



(13-4) 
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where j is any integer. 

(b) The constant zeta and omega loci in non-parametric form 
as developed in Section 5: 



f ) - 0 (13-5) 

W(a k , Id ) = 0 (13-6) 



(2) Other Specifications 

Consider the rational network function (13-7) which may be 
the transfer function of some network: 



z 

k=0 




t 

Lr = n 




(13-7) 



where the coefficients 
parameters q^: 






and b, 
k 



are different 



functions of the adjustable 



a k = f k* q l* q 2* ,,, * q i* (!3-8) 

b k = ®k^ q l’ q 2* * * * q i^ ( 13 ‘ 9 ) 

Specifications on the behavior of F(s) must be written as algebraic equations 
in the parameters q^ in order to apply this design technique. This is 
straightforward for the common types of specifications. Sections 6 and 7 
show how this is accomplished for specifications of Q and bandwidth. 

(3) Solution for the Parameter Values which Meet Specifications 
If r poles and zeros of the network function are specified, 
then Equations (13-3) through (13-6) may be employed as necessary to gener- 
ate r independent equations in the adjustable parameters q^. These r equa- 
tions taken together with perhaps s other equations arising from other 
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specifications form a set of (s+r) equations which must simultaneously be 
satisfied. It is presumed that the number of adjustable parameters 
is equal to or exceeds the number (s+r) of specifications. In order to 
obtain discrete solutions for the parameter values, the number of specifi- 
cations must be made equal to the number of adjustable parameters. This 
may be achieved by fixing as many of the adjustable parameters as required 
to obtain this equality. Once the number of adjustable parameters equals 
the number (s+r) of specifications, the set of (s+r) algebraic equations 
may be simultaneously solved. Simple cases may be solved by hand, however 
the equations are generally complicated and require computer solution. 
Several iterative procedures are available to achieve this, and are usual- 
ly some type of search routine which seeks parameter values which will force 
all equations to zero simultaneously. 

In general the set of (s+r) equations is non-linear and its solution 
will result in several sets of parameter values which mathematically satisfy 
all equations simultaneously. Some of these solutions will be rejected be- 
cause they are unrealizable. It should be kept in mind that the mathemati- 
cal solution may not always be physically acceptable even though realizable. 
This situation may occur, for example, when specifying a maximum in the 
frequency response characteristic by selecting parameter values which will 
cause the derivative of the magnitude function to be zero at some specified 
frequency. The mathematics may indicate a solution which results in a 
minimum rather than the desired maximum. In any event solutions should be 
checked to ensure that they meet all of the specifications. If a solution 
exists which is both realizable and physically meets the specifications, 
then it will be found among the solutions to the (s+r) equations. 
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d. Example 



Consider the case of the symmetrical parallel tee network 
which was investigated in Section 9. The zeros of the transfer function of 
this network are given by the roots of: 

of/3 p 3 +2(4p 2 +2/3p+l = 0 

Since this equation has two variable parameters, two of its three roots may 
be specified. A preliminary analysis has indicated that complex roots are 
possible with real and positive parameter values. Assume that the para- 
meter values are desired which will cause a pair of complex roots to have 
zeta equal to 0.2 and omega equal to 1.0. The basic parameter plane equa- 
tions for complex roots are chosen with j = 0 and j = 1. These become: 

0.84 © </3 -1.2/3 =0 
1.2 /S = 0.84 

Solution of these equations produce the desired parameter values: 

^ = 1.44 
jS = 0.70 

By dividing the quadratic factor associated with the specified complex roots 
into the cubic equation (13-10), the real root is determined to be located 
at p = -1. 

' e. Example 

Consider again the parallel-tee network of the previous example. 
Instead of completely specifying the complex zeros, let only the omega be 
specified, and instead of specifying zeta let the value of the real root be 
specified. Assume that the omega of the complex zeros is specified as 1.2 
and that the real zero is to occur at s 3 -0.5. In order to obtain the 
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desired parameter values which will result in these roots, the equations of 
the real roots and the constant omega loci are employed: 



Z \ <rk = 0 

k=0 

W ( a k’ = 0 

From Section 9, and with 00 - 1.2 and - -0.5, these equations are: 

-0.125 °</9 + 0.50/3 - 1.0/d + 1.0 = 0 (13-11) 

2.99 * 2 2 - 4.14 offi 2 + 2.88/3 - 1.0 = 0 (13-12) 

Equation (13-11) is solved for o( as a function of : 



o( = 0.50 /S - 1.0 

-0.125/3 



(13-13) 



Equation (13-13) is substituted into Equation (13-12) to form the quadratic 
in /3 : 



64.33/3 2 -221.34/6 + 190.10 = 0 



(13-14) 



Equation (13-14) is solved for two values of and then Equation (13-13) 
is used to obtain the corresponding values of °( . The results are: 



1 


= 0.85 


(3 y = 1.65 


CNI 

> 


= 0.46 


/* 2 " U79 



Upon checking of these results ^ and /$ ^ are found to produce the desired 
zeros. ^ and (3 ^ are rejected since they produce three real zeros, one 
of which has the desired sigma, rather than one real root and a pair of com- 
plex conjugate roots as desired. This situation was anticipated as discuss- 
ed in the introduction to this section. 
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f . Example 

In the previous two examples only zeros of the transfer func- 
tion were specified. In this example a pole is specified and a constraint 
is placed on the zeros. The symmetrical parallel tee network discussed in 
Section 9 is used as a vehicle. The transfer function of this network is: 



E o (p) 



c>V3p 3 + 2gp 2 +2/3p+i 



(13-15) 



E (p) *75 p + (2 / 2 > + 2 */3 + <* )pi(2 /3 + o( + 2)pfl 

where p ■ RCs 

Assume that a pole is required to be placed at p = -0.30 and that a pair of 
complex conjugate zeros is to have a damping ratio ^ = 0.4999. The equa- 
tion of the constant zeta loci for the numerator polynomial which was obtain- 
ed in Section 9 is used to express the constraint on zeta: 



((8.00 ©< -16.00) /3 +8.00)/3 = 0 



(13-16) 



The equation for the constant sigma curves for the denominator polynomial, 
also obtained in Section 9, is used to express the constraint on the pole: 



0.153*y3 -0.420/3 -0.050o< +0.400 = 0 

Equation (13-1-6) is solved for as a function of £ : 

o*' = 2.00 - (1.00//$ ) 



(13-17) 



(13-18) 



When Equation (13-18) is substituted into Equation (13-17) the following 
quadratic results: 



+ 1.52/S - 1.84 = 0 



(13-19) 



Equations (13-18) and (13-19) give the following pair of parameter values 
which mathematically satisfy the constraints placed on the pole and zeros 
of the transfer function: 



* l = 0.74 
2 = 2.43 



/ 3 . 



/*2 = - 



0.80 



= -2.31 
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Upon checking, and ^3 ^ are found to satisfy the specified constraints. 

0( 2 and f$ 2 are rejected as unrealizable since /S ^ is negative. 

g. Example 

Consider the loaded and null-adjusted symmetrical parallel 
tee network discussed. Presume that one wishes to select a value of the 
source resistance &( , and the load resistance (3 , which will cause a 
pole to be located at p ■ 0" * -0.3 and which will result in a value of 
Q :^0.34. Referring to Section 10, Equation (10-2), it is seen that the 
specification on the real pole requires that: 

-0.84 4(3 - 0.11 -0.11 /3 + 1.40 * 0 (13-20) 

Equation (10-10) of the same section gives the equation which must be satis- 
fied to meet the specification on Q: 

1.8496 * 2 (S 2 - 0.3008 2 + 0.8496 £ 2 - 0.3008 o < 2 /S - 4.4512 - 

0.1504(9 + 0.8496-< 2 - 0.1504y + 0.4624 = 0 (13-21) 

Equations (13-20) and (13-21), when simultaneously solved by digital computer 
produce the solution: 

* = 1.95 

(3 = 0.677 

These values are physically realizable and meet the specifications. 

h. Conclusions 

The parameter plane technique presented in this section is 
the only design technique known to the author which allows parameter values 
to be determined which will cause both pole-zero and performance specifica- 
tions to be simultaneously satisfied. The types of specifications which 
may be made are limited only by the requirement that they be expressible as 
algebraic equations in the adjustable parameters. For common specifica- 
tions this is readily achieved. 
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14. Conclusions and Recommendations for Further Work 



a. Form of the Coefficients 

The work which has been presented in Section 4 shows how parameter 
plane equations have been developed for polynomials with coefficients of 
the form: 

\ - V* + + Vi 3 + £ k f 14 - 1 ) 

This extension of the theory permits a wide class of network functions to 

be analyzed. There remains, however, cases where higher order terms in 

2 2 2 

and appear in the coefficients. These include terms like: cx /S , «< ft , 

2 

cxft , etc. Work currently being performed by Lt. A. R. Miller, USN at 
the U. S. Naval Postgraduate School has shown that it is possible to solve 
the parameter plane equations for coefficients of the form: 

\- b k ° <2 + 'kP 2 + Vf 3 + v< +£ kP + 3 k <14 - 2) 

Preliminary work done by Lt. R. E. Hudson, USN, also at the U. S. Naval 
Postgraduate School, has shown that it is theoretically possible to solve 
the parameter plane equations for coefficients of the form: 

•k ■ v* 2 ** 2 + c k* 2 0 + “k 2 + e k°< 2 + w 2 + w 8 + \« + ■>ki 3+ c k 

(14-3) 

There remains however much to be done, particularily in interpreting the 
curves which result from solution of the parameter plane equations associated 
with (14-2) and (14-3). 

It is, of course, desirable to allow even higher order combinations of 
o( and /S to appear in the polynomial coefficients than allowed by (14-3). 

At the current state of development it is not clear just how to achieve this. 
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b . Non-parametric Equations of the Constant Zeta and Omega Loci 

The derivation of the non-parametric equations of the constant 
zeta and omega loci presented in Section 5 permits these loci to be more 
readily constructed than was previously possible. More important, however, 
it forms the bridge between parameter plane theory and the theory of plane 
algebraic curves. To date this relationship has not been seriously investi- 
gated. It is the author 1 s opinion that such investigation will lead to new 
and valuable information about the behavior of the loci. Such information 
may well lead to the development of sketching techniques which would con- 
siderably simplify the construction of the loci. Additionally it might 
permit new insight into the behavior of the zeta equals zero locus which is 
the stability limit for the poles of active networks, 
c. Constant Bandwidth Curves 

Simultaneous correlation between frequency response character- 
istics and root locations upon the parameter plane is now possible as a 
result of the development of the constant bandwidth curves which were develop- 
ed in Section 6. Families of such curves permit the frequency response of a 
network to be sketched for any choice of parameter values within the range 
for which the curves are drawn. 

The equations developed in Section 6 apply to coefficients in the trans- 
fer function of the form (14-1). When parameter plane techniques are 
developed for coefficients of the form (14-3) it will be necessary to extend 
the work of Section 6 to include this coefficient form. Preliminary work 
done by the author has shown that this is possible by following the same 
procedure as employed in Section 6. 
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d. 



Constant Q Loci 



As defined by Morris (4), Q is a mathematical parameter of a 
second order characteristic equation. Thus defined, Q may be employed to 
characterize RC networks for which the conventional definitions fail to 
provide useable correlation between the network characteristics and the 
value of Q. 

Loci of constant Q, developed in Section 7, permit the value of Q 
associated with any choice of network parameters to be read directly from 
the parameter plane. Alternately one may select parameter values from the 
parameter plane which will achieve a desired value of Q within the capabili- 
ties of the network. 

e . Algebraic Design by Parameter Plane Methods 

The algebraic design techniques presented in Section 13 per- 
mit m network parameter values to be determined which will guarantee that m 
specifications upon the network are met. These specifications can be either 
pole-zero locations or other specifications which can be written as algebraic 
equations in the m adjustable parameters, or both. As pointed out in 
Section 13, it is necessary to have the same number of adjustable parameters 
as specifications in order to obtain discrete sets of parameter values as 
solutions to the problem. Unless these specifications are realistic ones 
in terms of the network capabilities, then it may be the case that no 
physically realizable solutions will be produced by the design technique. 

It is more likely that a realizable solution would be found if the system 
had more adjustable parameters than specifications and if the specifications 
were stated as inequalities rather than strict equalities. For example, one 
might specify Q greater than 0.34 rather than strictly equal to 0.34 provid- 
ing that this is acceptable. This procedure would lead to a set of equations 
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which would contain some strict equalities and other inequalities and hence 
other methods than those of Section 13 would have to be applied in order to 
obtain parameter values which would meet the specifications. The author 
suggests that non-linear programming techniques be investigated in connec- 
tion with the solution of such a set of equations, 
f . Summary 

The basic parameter plane techniques, and the extensions to 
them presented here, provide a powerful tool for the analysis and design of 
variable parameter, active and passive, electric networks. These techniques 
are capable of providing information about the poles and zeros and sensitivity 
of variable parameter networks which is not possible by any other method 
currently known. It is, for two parameter problems, all and more that 
Evan's root locus technique is for single parameter problems. Like the 
root locus technique it is applicable to a much wider sphere of problems 
than that of electric networks. Subject to the restrictions on the form of 
the coefficients previously discussed, parameter plane techniques may be 
used to investigate the singularities of any variable coefficient rational 
function in the complex variable s. 
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APPENDIX I 



Inductive proof that: 

u i< j ? )u j ( r > - - u i. ,< > < i - i > 

for all integers i and j. 

The inductive proof that (I - 1) is true consists of the following steps: 

1. Show that (1-1) is true for i = 2 and i = 3. 

2. Assume that (1-1) is true for i = m-1 and i = m. 

3. Show that (1-1) is true for i = m + 1, given that it is true 
for i = m-1 and i = m-2. 

Step la: 

1 = 2 

U 2 ( f )U j-l<? > - U i<? )Uj<? ) = _U 2-j ( f ) d-2) 

but 

u 2 ( f ) = 2 ^ 

V ? > - 1 

Therefore (1-2) may be written: 

2 ^ U j-1 ( 7 ) - U j ( f ) - - u 2 .j< T > (I_3) 

but U_ k ( £ ) - -U k ( f ) 

Hence (1-3) becomes: 

j 7 u j-i < P‘V* ) - u j- 2 < ? > (i - 4> 

(1-4) is recognized as the recursion relation which defines the 
Chebyshev functions. Thus (1-2) is true for all j. 

Step lb: 
i = 3 

B 3< f > U j-1< ? > ' B 2 < 5 >V T > - '“a- J< ? > (I - 5) 
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but : 



U 3< ? > ' 4 f - 1 
Therefore (1-5) may be written: 

(4^ 2 -l)U j :L (^ ) - 2 ^ U.( f ) = U._ 3 ( f ) (1-6) 

Making use of the recursion relation for the Chebyshev functions, (1-4), 
it is possible to write (1-6); 

(4^ 2 -l)U j ;L ( £ ) - 4^ 2 U._ 1 ( ^ ) + 2^U._ 2 ( ^ ) = U._ 3 ( J ) (1-7) 

Expanding (1-7) and cancelling where possible^ produces the recursion relation 
for the Chebyshev functions, (1-4). 

Therefore (1-5) is true. 

Step 2: 

Assume that the following equations are tjjj?e: 

U ~l<? ) Vl < y > • V 2 <?> U j<5 ) ■ > < I ' 8) 

V ? > Vi< 7 > - Vi ( ? )u j ( ? > ■ -Vj< ? ) < M > 

Step 3: 

v>n< ? > ■ w y > (i - io> 

Applying the recursion relation for the Chebyshev functions, (1-4), to the 
left hand side of (1-10) and rewriting, one obtains: 

-Vj + i<P (I - U) 

Substituting (1-8) and (1-9) into (1-11) one obtains: 



-u 



m-j+l( £) 



( 1 - 12 ) 



Equation (1-12) may be written: 
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< 2 5 U »<f > ‘ U *-l < ? ),U j-l ( P - [2 f U n,.l<P-V2<f> 1 U j<?> ' 

(1-13) 

By use of the recursion relation for the Chebyshev functions, (1-4), Equa- 
tion (1-13) may be written: 

• V?> u j<5> ' (I - M> 

Thus (1-1) is true for i = m+1, given that it is true for i = ra-1 and i = 
m-2. This concludes the proof that (1-1) is true for all i and j integer. 



97 



APPENDIX II 



Proof that: 



t. ( ‘ 1 ) \" ku k-j ( f ) ’ 0 



k=0 



(II-D 



for all integers j and k 



Theorem: 



If: ^ (-l) k a, AJ k U, ( X ) = 0 

k=0 



and if: 



t <- i > V\-i ( ? > = 0 

k=0 



(II-2) 



(II-3) 



Then: 



irs 



(-l) k a^to k U^ j( ^ ) = 0 for all integei 
k=0 j and k (II -4) 



Proof: 



Subtract (II-2) multiplied by ^ ) from (II-3) multiplied by 

U.( ^ ). This gives: 



n 

Z <- 1 >'V k • ’WWP 1 -° 



(II-5) 



k=0 



Employing the results obtained in Appendix I, namely (1-1), Equation (II-5) 



may be written: 

r\ 



J. -0 



(II-6) 



k=0 

Multiplying (II-6) by -1 concludes the proof of (II-4). 
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APPENDIX III 



Table of the Chebyshev Functions 

The Chebyshev Functions T^( ^ ) and U^( ^ ) are defined by the following 



recursion relations: 



W? > - 2 + > -° 

v ? > - 1 

V 5 > ' 

W 3 > -W 5 > + \-i< 3 > ’ 0 

V 3 > * 0 

V 3 ) ■ 1 

The first few of these functions are tabulated below: 



k 




\< 3 > 



0 

1 

2 

3 

4 



1 



3 

2 5 2-1 

45 3 -3J 
8 £ 4 -8 J 2 +l 



0 



2 ? 

4? 2 -l 

8 3 3 - 4 3 
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